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ABSTRACT 


The  idea  of  survivability  introduced  as  a  network  parameter  has  led 
to  so  many  investigations.   Several  measures  of  survivability  has  been 
studied.   The  number  of  links  and/or  stations  needed  to  be  damaged  to 
disrupt  the  system  is  the  survivability  criterion  adapted  in  this  study. 

The  development  of  analysis  procedures  for  directed,  undirected,  or 
mixed  networks  based  on  the  above  criterion  and  use  of  the  concepts  in 
network  flow  and  graph  theory  are  treated  in  detail  including  computer 
program  implementation  of  the  algorithms.   Finally  a  practical  design 
algorithm  for  minimum-cost  survivable  network  with  respect  to  branch  dis- 
connection using  a  heuristic  approach  and  analysis  techniques  is  described. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 6 

II .  NETWORK  FUNDAMENTALS 8 

2. 1  NETWORK  TOPOLOGY 8 

2.1.1  Directed,  Undirected,  and  Mixed  Networks 8 

2.1.2  Paths  in  a  Network  Graph „ ..  .  .  9 

2.1.3  Cutsets  and  Cuts 9 

2.1.4  Network  Tree 11 

2.2  FLOWS  IN  DETERMINISTIC  NETWORKS 11 

2.2.1  Branch  Flow  and  Branch  Capacity 11 

2.2.2  Flow  Pattern 12 

2.2.3  Residual  Capacity  and  Augmentation  Path 12 

2.2.4  Maximizing  Network  Flow. 13 

2.3  MATRICES  ASSOCIATED  WITH  SURVIVABLE  NETWORKS 14 

2.3-1  Incidence  Matrix...... ..14 

2.3«2  Connection  Matrix 14 

2.3.3  Branch  Capacity  and  Terminal  Capacity  Matrices. ...  14 

2.3.4  Connectivity  and  Redundancy  Matrices .....15 

2.4  DISCONNECTING  SETS  OF  A  NETWORK 15 

III .  ANALYSIS  OF  DETERMINISTIC  SURVIVABLE  NETWORKS 17 

3.1   SURVIVABILITY  CRITERIA 0....0 17 

3.  2  ANALYSIS  WITH  RESPECT  TO  BRANCH  DISCONNECTION 18 

3.2.1  s-t  Branch  Connectivity , 18 

3.2.2  Multi terminal  Branch  Connectivity 22 

3.3  ANALYSIS  WITH  RESPECT  TO  VERTEX  DISCONNECTION 32 

3.3.1   s-t  Vertex  Connectivity .32 

3 


3.3.2     Multiterminal  Vertex  Connectivity..... 36 

IV .  DESIGN  OF  LOW-COST  SURVIVABLE  NETWORK 39 

4.1  INTRODUCTION 39 

4.2  DESIGN  WITH  BRANCH  DISCONNECTION  CRITERION 39 

4.2.1  The  Starting  Routine 40 

4.2.2  Optimizing  Routine . .  .41 

4.2.3  Feasibility  Test 44 

V.  CONCLUSION ' 50 

APPENDIX  A     FLOWCHARTS . 51 

A-l      s-t  BRANCH  CONNECTIVITY 51 

A-2     MULTITERMINAL  BRANCH  CONNECTIVITY 

(Undirected  Network) 52 

A-3     MULTITERMINAL  BRANCH  CONNECTIVITY 

( Directed  Network) 5^ 

A-4     s-t  VERTEX  CONNECTIVITY 56 

A-5     UNIFORM  R-CONNECTIVITY 57 

APPENDIX  B     COMPUTER  PROGRAMS 58 

B-l     MAXCON-1 58 

B-2     MAXCON-1  A 64 

B-3     MAXCON-2.. 67 

B-4     MAXC0N-2A 76 

B-5     MAXCON-3 78 

B-6     MAXC0N-3A 84 

B-7     MAXCON-4 87 

BIBLIOGRAPHY 91 

INITIAL  DISTRIBUTION  LIST 93 

FORM  DD  1473 0 94 


ACKNOWLEDGEMENT 

I  wish  to  express  my  gratitude  to  my  thesis  advisor,  Doctor  Charles 
H.  Rothauge  of  the  Department  of  Electrical  Engineering,  Naval  Postgrau- 
ate  School,  for  helping  me  strengthen  my  background  in  network  flows,  and 
for  his  continued  guidance  and  encouragement  during  all  the  phases  of  my 
work. 


-X 


I.   INTRODUCTION 

Network  complexes  such  as  railroads,  highways,  communication  circuits, 
telephone  systems,  and  many  others  which  are  interconnections  of  links  and 
stations,  and  where  there  exist  commodity  flow,  have  been  modelled  by  lin- 
ear graphs.   The  vertices  of  the  graph  symbolize  stations  or  terminals 
while  the  branches  denote  lines,  trunks,  or  circuits.   The  graphical  struc- 
ture of  such  systems  allows  the  use  of  topological  concepts  in  the  develop- 
ment of  the  theory  of  network. 

In  1961  the  idea  of  survivability  as  a  graph  parameter  was  introduced. 
This  concept  was  later  adapted  as  a  network  parameter.   Survivability  of 
a  network  is  defined  as  its  capability  to  function  as  a  system  after  dam- 
age induced  by  enemy  attack.   Survivability  and  vulnerability  are  used 
interchangeably  by  many  authors.   The  term  network  reliability  is  more  ap- 
propriately used  when  referring  to  damage  caused  by  natural  disturbances. 

The  study  of  network  survivability  employs  topological  methods  and  the 
concept  of  network  flow.  There  are  two  distinct  areas,  namely:   determin- 
istic survivability  and  probabilistic  survivability.   A  network  is  con- 
sidered probabilistic  when  the  parameters  involved  like  commodity  flow  and 
weights  assigned  to  links  and  stations  are  statistical.   When  there  is 
certainty  in  the  existence  of  these  parameters,  the  network  becomes  de- 
terministic. 

Up  to  present,  network  survivability  has  been  the  subject  of  research 
by  many  graph  and  network  theoreticians.  Recent  works  were  focused  on  the 
mathematical  formulation  of  survivability  criteria,  the  development  of 
analysis  methods,  said  the  synthesis  of  optimally  invulnerable  networks. 
Numerous  results  have  already  been  published. 


This  thesis  is  based  on  the  paper  by  Frank  and  Frisch  [lO^j  which 
summarizes  the  most  significant  results  of  the  researches  in  this  area. 
The  treatment  of  deterministic  networks  is  restated  in  more  detail  and  ex- 
tended by  adapting  other  results  in  network  flow  theory  with  some  modifi- 
cations to  provide  additional  analysis  procedures.   To  provide  working 
tools  for  analysis  and  synthesis,  the  existing  algorithms  and  those  that 
were  modified  were  coded  for  the  computers.   Further  simplifications  were 
introduced  when  encoding  these  algorithms  into  the  programs  which  appear 
in  the  appendices  attached  herewith. 

In  Chapter  II  of  this  report,  a  summary  of  the  important  elementary 
principles  of  network  theory  essential  in  the  study  of  network  surviva- 
bility is  presented.   Chapter  III  covers  the  analysis  procedures  based  on 
some  survivability  criteria.   The  principles  from  which  they  were  derived 
are  also  discussed.   Investigation  of  directed,  undirected,  and  mixed  net- 
works with  respect  to  branch  disconnection  and  vertex  disconnection  are 
separately  treated. 

These  analysis  tools  may  be  employed  in  the  design  of  optimally  sur- 
vivable  networks.  The  object  of  Chapter  IV  is  design  optimization  by  an 
heuristic  approach  using  analysis  techniques. 


II.   NETWORK  FUNDAMENTALS 

2.1  NETWORK  TOPOLOGY 

One  of  the  many  applications  of  linear  graphs  to  engineering  problems 

is  the  topological  representation  of  networks.  Networks  may  be  categorized 

as  "electrical  networks",  "switching  networks",  and  "flow  networks".   The 

latter  which  is  often  called  "communication  net"  has  associated  with  its 

branches  some  information  or  commodity  flow.   Symbols,  terminology,  and  a 

minimum  of  the  fundamental  concepts  necessary  in  the  study  of  this  type  of 

network  is  discussed  in  this  chapter. 

2.1.1   Directed,  Undirected,  and  Mixed  Networks 

A  directed  network  is  represented  by  a  directed  linear  graph  G=(V,r), 

consisting  of  a  set  of  elements  called  vertices  denoted  by  the  symbol  V 

and  a  set  of  ordered  pairs  of  vertices  called  branches  denoted  by  P  .  The 

vertex  v.  is  called  the  initial  vertex  and  the  vertex  v  .  the  terminal  ver- 
i  J 

tex;  vertices  v.  and  v.  are  said  to  be  adjacent. 

A  directed  branch  is  indicated  by  the  ordered  pair  (i,j).   Each 
branch  carries  an  orientation  and  the  information  flow  is  only  in  the  spec- 
ified direction. 

An  undirected  network  is  represented  by  an  undirected  graph  G-(V ,V) , 
consisting  of  a  set  of  vertices  V,  and  a  set  of  unordered  pairs  of  these 
vertices  denoted  by  p  .   An  undirected  branch  is  indicated  by  [_i»jj  or 
LJ,i^|.   The  flow  of  information  is  bidirectional. 

A  mixed  network  is  symbolized  by  a  graph  that  has  both  directed  and 
undirected  branches.  A  non-oriented  branch  may  be  represented  graphically 
by  two  parallel  branches  oppositely  directed. 


The  network  graph  may  have  more  than  one  branch  between  the  same 
pair  of  vertices. 

2.1.2  Paths  In  a  Network  Graph 

A  directed  s-t  path  tt.  is  a  sequence  of  distinct  vertices  and 
branches  connecting  two  arbitrary  vertices,  v  and  v,  .   It  has  the  form 

v  (s'V'  V  ^l'3^-  vi  ■  •••  '  v  (ik,t)l  Vt 

1  <c  K 

An  undirected  s-t  path  n.  is  a  sequence  of  distinct  vertices  and 
branches  of  the  form 

V  [>'V1'V  E1!'^'  vi  '  —  •  Vi  '  [ik't]'  Vt 

1  *c  K 

Two  directed  or  undirected  paths  tt.  and  tt  .  are  branch  disjoint  if 

they  have  no  branches  in  common  and  vertex  disjoint  if  the  only  vertices 

they  have  in  common  are  v  and  v, . 

s      x» 

In  a  sequence  of  vertices  and  branches  describing  a  directed  path, 
if  (i, j)  is  in  the  path,  then  (i,j)  is  a  forward  branch  with  respect  to 
the  path.   A  branch  (j,i)  in  the  path  is  a  backward  branch  with  respect 
to  the  path. 

The  network  graph  G  is  said  to  be  connected  if  for  any  pair  of  ver- 
tices v.  and  v.  there  is  a  path  tt.  .  in  G.   Otherwise  it  is  disconnected. 
The  resulting  subgraphs  are  components. 

2.1.3  Cutsets  and  Cuts 

An  undirected  branch  cutset  of  an  undirected,  connected  graph  is 
the  minimal  set  of  branches,  the  removal  of  which  yields  a  graph  of  two 
or  more  components. 

A  directed  branch  cutset  is  a  minimal  set  of  branches  of  a  directed 
graph,  the  removal  of  which  disconnects  all  directed  paths  from  at  least 
one  vertex  of  G  to  at  least  another  vertex  of  G. 


An  undirected  vertex  cutset  of  an  undirected,    connected  graph  is 
the  minimal   set  of  vertices  whose  removal  from  G  disconnects  the  graph. 

A  directed  vertex  cutset  is  a  set  of  vertices  of  a  directed  graph 
whose  removal  from  G  destroys  all   directed  paths  from  at  least  one  of  the 
remaining  vertex  of  G  to  at  least  another  remaining  vertex. 

Cutsets  with  respect  to  a  vertex  pair  v     and  v,    are  referred  to  as 
s-t  cutsets. 

A  cut   (X,X)    separating  v     and  v,    is  a  set  of  branches  (i,j)    such 

S  °t> 


that  v     €  X  and  v     G   X.      In  Fig.    2.1,    if  X 

S  w 


J v    ,v    ,v    I     and  X  -     ( v    , 
vk»v+  1     i    the  set  of  "branches   (X,X)   ■      {  (i,^),    (1,3),    (2,3)  >     is  a  directed 


s-t  cut.      Branch  (2,^4-)    is  oppositely  directed  with  respect  to  the  vertices 
s  and  t  therefore  it  does  not  belong  to  the  set. 


Fig.    2.1 
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2.1.4     Network  Tree 

A  spanning  tree  T  =  (V,U)    of  a  connected  network  graph  G  =   (V,  P   ) 
is  a  connected  subgraph  of  G  containing  all  of  the  vertices  of  G  hut  no 
closed  path. 

If  n  is  the  number  of  vertices  in  G,    a  tree  contain  n-1  branches. 

Every  pair  of  vertices  v.    and  v.   of  graph  T  is  joined  by  one  and 
only  one  i-j  path.     Therefore  it  becomes  disconnected  when  one  branch  is 
removed. 

2.2     FLOWS  IN  DETERMINISTIC  NETWORK 
2.2.1     Branch  Flow  and  Branch  Capacity 

Aside  from  the  topological   structure  of  a  communication  net,   there 
is  assigned  to  each  branch  a  finite  weight  f(i,j)    or  f[i,j^  called  the 
branch  flow  from  v.    to  v..      Also  a  weight  c(i,j)   or  c[i,j]  is  assigned  to 
each  branch  which  is  the  maximum  amount  of  commodity  that  can  reach  v. 

from  v.    per  unit  time,    called  the  branch  capacity. 
For  an  undirected  network, 

c[i,j]  =  -c[j,i],  (2.2.1.1) 

f[i,j]  =  -f[j,i],  (2.2.1.2) 

f[s,t]  =  I  f[i,j],  (2.2.1.3) 

where  v     and  v,    are  any  chosen  source  and  sink  vertices  respectively.      For 

S  w 

a  directed  case,  equations  (2.2.1.1)  and  (2.2.1.2)  are  not  necessarily 
true. 

An  undirected  branch  can  be  represented  by  two  oppositely  directed 
branches  between  v.  and  v.  with  capacities 

c(i,j)  -  c(j,i)  =  c[i,j],  (2.2.1.4) 
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2.2.2  Flow  Pattern 

The  set  of  flows  associated  with  the  branches  in  G  is  called  a  flow 
pattern  denoted  "by  F,   A  flow  pattern  is  said  to  be  feasible  if  it  satis- 
fies Equations  (2.2.2.1)  and  (2.2.2.2)  for  some  nonnegative  constant  f 
which  is  called  the  value  of  F. 

For  all  v.  €  V  of  a  directed  network, 
l 


' f  .   if  i=s 
s,t 


f(i,V)  -  f(V,i)  =, 


0     if  i^s,t  (2.2.2.1) 


f  .   if  i=t 
V.  s,t 

0  <  f(i,j)  <  c(i,j)  for  all  i  and  j         (2.2.2.2) 

This  means  that  the  net  flow  out  of  the  source  v  is  f   ,  and  the 

s     s,t 

flow  out  of  the  sink  v,  is  -f   ,  ,  whereas  the  net  flow  out  of  an  intermed- 

t      s,t' 

iate  vertex  is  zero.   This  latter  relation  is  also  true  for  an  undirected 

network. 

A  branch  (i,  j)  is  saturated  if  f(i, j)  =  c(i, j)  and  unsaturated 

otherwise.   It  is  said  to  be  empty  if  f(i,j)  =  0. 

2.2.3  Residual  Capacity  and  Augmentation  Path 

The  residual  capacity  r(rr  ,)  is  defined  both  for  directed  and  un- 

s  i  x> 

directed  network  as 

r(TTg  t)   =  Min       Min    [f(i,j)J     ,   Min    [c(a,b)    -  f(a,b)] 
'  L  backward  forward  J 

branches  branches 

(i,j)   €  ns,t  (a,b)   €  ns,t  (2.2.3.1) 

A  path  n  .  is  saturated  if  r(n  ,  )  =  0  and  unsaturated  otherwise. 
*  s,t  v  s.t-^ 

Additional  flow  can  be  sent  through  an  unsaturated  path  obtaining  a  new 
flow  pattern.   To  achieve  this,  add  r(n   ,)  to  all  forward  flows  and 
subtract  r(rr  ,)  from  all  backward  flows  along  the  path. 
If  a  path  between  two  vertices  has 
f(i,j)  <  c(i,j) 
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and  f(j.i)  >  0, 

then  we  call  this  a  flow  augmentation  path.   The  flow  along  this  path  can 

he  increased  according  to  the  previously  stated  rule. 

2.2.4  Maximizing  Network  Flow 

One  of  the  most  important  and  "basic  theorems  in  network  flows  is 

the  "Max- flow  Min-cut  theorem"  formulated  by  Ford  and  Fulkerson  [_7~]»      The 

problem  of  obtaining  the  maximum  obtainable  flow  value  f    through  a  net 

s  j  x 

with  reference  to  two  vertices  v     and  v,    is  solved  using  this  theorem. 

s  t 

Max-flow  Min-cut  Theorem:      The  maximal  flow  value  obtainable  in  a 

network  G  form  v     to  v,    is  the  minimum  of  the  value  of  a  cutset  taken  over 
s  t 

all   cutsets  separating  v     and  v , . 

To  state  the  theorem  mathematically,    let: 

f       .    be  the  value  of  the  minimum  s-t  cut 

c(X. ,X. )    be  the  capacity  of  a  cut   (X. ,X. ) 


{<W..t} 


be  the  set  of  all  s-t  cuts  in  the  net 


then: 


Max  [f   ]=  Min  [c(X  ,X  )   ]    =T   ,     (2.2.4.1) 

The  value  v  of  the  cut  (X,X)  is  defined  as 

v  =  f(X,X)  -  f(X,X)  <  c(X,X)  (2.2.4.2) 

The  Ford  and  Fulkerson  Labeling  and  Augmentation  algorithm  [y~\  based 

on  the  above  theorem  provides  a  method  of  increasing  f      .    systematically 

S  i  *c 

from  an  arbitrary  feasible  flow  pattern  with  value  less  than  f      ,  and  ter- 

S  I  v 

minates  when  f   ,  =  f      , .   The  resulting  flow  pattern  is  maximal. 

The  algorithm  will  be  stated  in  a  modified  form  in  Chapter  III. 
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2.3     MATRICES  ASSOCIATED  WITH  SURVIVABLE  NETWORKS 

Network  graphs  may  "be  structurally  drawn  or  be  described  by  matrices. 
The  latter  representation  is  useful  for  computer  calculations.      This  sec- 
tion enumerates  several  matrices  used  in  the  study  of  survivable  networks. 
2.3.1     The  Incidence  Matrix 

Let  G  =  (V,  P  )  be  the  directed  graph  of  the  network  consisting  of 
n  vertices  and  m  branches.  Let  A  =  fa.  .1  be  an  n  x  m  matrix  whose  i-ith 
entry  is  a.     ..      Then  A  is  the  incidence  matrix  of  G  if 

+1     if  a  branch  is  directed  away  from  v. 


a.  .  =  < 
if  J 


I 


-1  if  a  branch  is  directed  towards  v. 

1 

0  otherwise 


If  G  is  an  undirected  graph,  A  =  [a.  .J  is  the  incidence  matrix  if 

1  >  J 


1     if  branch  is  incident  at  v. 

1 

a.     .  =    1 

1    0     otherwise 
2.3.2     Connection  Matrix 

An  alternate  representation  of  a  non-oriented  network  is  the  con- 
nection matrix  K  =  [k.     .J,      K  is  an  n  x  n  symmetric  matrix  whose  i-jth 

1 » 3 

entry  is  k.     .  where 

1      if  [i,j]  €  P 


k.     . 
i»3 


0     otherwise 


For  a  directed  network,    a  -1  is  entered  for  an  oppositely  directed  branch. 

2.3.3     Branch- capacity  and  Terminal   Capacity  Matrices 

The  other  way  of  representing  a  communication  net  of  n  vertices  is 

by  using  a  symmetric  matrix  C  =  \_c.     .J  called  the  branch- capacity  matrix 

1 » J 

where  c.     .  =  c .    .    is  the  sum  of  the  capacities  of  the  branches  connected 
between  v.    and  v..     For  an  oriented  network,    the  C-matrix  is  asymmetrical. 

14 


For  any  two  vertices  in  G,  there  exist  a  finite,  maximum  obtainable 
capacity  of  transferring  information  from  one  vertex  to  the  other  called 

the  terminal  capacity.   The  terminal  capacity  matrix  T  =  [t.     .j  represents 

1  *  3 

the  maximum  flow  that  can  exist  between  all  the  possible  vertex  pairs  of 
G.     The  T-matrix  is  also  asymmetrical  for  a  directed  network. 
2.3.^     Connectivity   and  Redundancy  Matrices 

The  connectivity  matrix  is  essentially  a  T-matrix  with  each  branch 
of  G  or  G  assigned  a  maximum  capacity  of  unity. 

An  n  x  n  matrix  R  =  [_t.     .1  is  a  redundancy  matrix  where  the  i-jth 
entry  r.     .is  the  number  of  branch/vertex  disjoint  paths  between  v.    and 
v..     It  is  symmetrical  for  a  non-oriented  net.      A  uniform  R-matrix  is  one 
which  has  identical   entries   except  the  main  diagonal. 

The  main  diagonals  of  these  matrices  have  infinity  as  entries  how- 
ever for  computer  calculation  zero  entries  are  used. 

2.k     DISCONNECTING  SETS  OF  A  NETWORK 

A  connected  network  represented  by  G  can  be  disconnected  by  deleting 
certain  branches  or  vertices. 

An  s-t  branch  cutset  is  a  set  of  branches  whose  deletion  breaks  all 

directed  paths  from  v     to  v , ;    likewise  an  s-t  vertex  cutset  is  a  set  of 

vertices  whose  deletion  breaks  all  directed  paths  from  v     to  v , .      A  mixed 

st 

s-t  cutset  is  composed  of  branches  and  vertices. 

Clearly  if  the  capacity  of  each  branch  or  vertex  is  set  to  unity, 
then  the  value  of  the  cutset  is  numerically   equal  to  the  number  of  com- 
ponents of  the  cutset. 


set; 


Let    Ci      ,    be  the  number  of  elements  in  the  smallest  s-t  vertex  cut- 
s,t 


f      ,    be  the  number  of   elements  in  the  smallest  s-t  branch  cutset; 
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C       ,    "be  the  number  of  elements  in  the  smallest  s-t  mixed  cutset; 
s,t  ' 

then  the  disconnecting  sets  of  the  connected  graph  are  as  follows: 

6)  =  Min  [U       ]  (2.4.1.1) 

S  j  X 


(2.4.1.2) 


s,t  €   V 

r  = 

Mln[TSjt] 

s,t  €    V 
s/t 

<r  = 

Min[Cs>t] 

s,t  €    V 
s/t 

(2.4.1.3) 


The  smallest  disconnecting  set  of  a  connected  graph  denoted  by  v  is 
defined  as 

v  -  Min  [  4J  ,  r  ,  <T    ]  (2.4.1.4) 

However  it  was  shown  in  \_2>~\   that 

v  =W=  Min  [  W   ,].  (2.4.1.5) 

S  f  X- 

s,t  6  v 
s^t 


16 


III.   ANALYSIS  OF  DETERMINISTIC  SURVIVABLE  NETWORKS 

3.1  SURVIVABILITY  CRITERIA 

A  significant  criterion  in  network  analysis  is  its  survivability  or 
vulnerability.  In  a  communication  net,  this  parameter  refers  to  the  de- 
gree in  which  the  system  remain  functional  after  an  enemy  attack  or  nat- 
ural disaster. 

No  general  survivable  criteria  has  yet  been  formulated  however  several 
measures  of  survivability  has  been  investigated  namely: 

1.  The  number  of  links  and/or  stations  that  must  be  destroyed  before 
communication  is  disrupted  {_^~]t 

2.  the  minimum  number  of  links  that  must  be  deleted  in  order  to  iso- 
late a  group  of  stations  from  communication  with  the  remainder  of  the  net 

W. 

3.  the  minimum  number  of  isolated  subnets  that  would  result  from  an 
attack  or  destruction.   This  is  termed  the  "independence  number"  of  the 
network  \_1~\, 

4.  fraction  of  stations  which  survive  an  attack  [2],  \_6~], 

5.  fraction  of  stations  that  can  be  reached  from  a  given  point  by 
a  path  of  no  more  than  r  links  after  an  attack  [8], 

The  first  three  criteria  apply  to  deterministic  network  while  the 
last  two  criteria  are  associated  with  probabilistic  net. 

Analysis  of  large  networks  based  on  the  first  criterion  is  treated 
in  this  thesis.  This  index  of  survivability  suggests  the  application  of 
the  concept  of  disconnecting  sets.  More  specifically,  the  parameter  (a) 
is  a  measure  of  vulnerability  to  station  destruction  and  f  is  a  measure 
to  link  destruction. 
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3.2  ANALYSIS  WITH  RESPECT  TO  BRANCH  DISCONNECTION 
3.2.1   s-t  Branch  Connectivity 

The  problem  discussed  here  is  the  finding  of  the  maximum  number 
of  branch  disjoint  s-t  paths  between  any  given  pair  of  vertices  v  and 

o 

v,    in  G.      Setting  the  branch  capacities  to  unity  and  applying  the  max- 

flow  min-cut  theorem,    the  maximal  flow  f      ,    from  v     to  v,    is  obtained 

s,t       s     t 

which  is  numerically  equal  to  the  smallest  branch  cutset  f        .   This 

Sj  t 

value  further  represents  the  maximum  branch  disjoint  paths  between  the 
given  vertex-pair. 

The  Ford-Fulkerson  algorithm  solves  a  general  max-flow  problem  be- 
tween two  vertices,  the  source  and  sink.  It  is  stated  here  in  simplified 
form  to  adapt  it  to  the  special  case  being  investigated. 

A  few  terms  are  defined  below  before  stating  the  algorithm: 

1.  A  vertex  is  labeled  by  an  ordered  pair  (i,+).  For  a  particu- 
lar vertex  v.,  the  index  i  denotes  the  preceeding  labeled  vertex  v. . 
The  second  term  in  the  label  is  plus  if  for  a  forward  branch,  f(i,j)  =  0 
and  a  minus  sign  if  for  a  backward  branch,  f(j,i)  >  0.   Otherwise  a  ver- 
tex is  unlabeled. 

2.  A  particular  vertex  v  is  scanned  if  its  adjacent  vertices  v 

are  labeled  from  it.  Otherwise  v  is  unscanned. 

x 

3.  A  vertex  is  considered  to  be  in  one  of  the  three  states,  namely: 

a.  unlabeled  and  unscanned 

b.  labeled  and  unscanned 

c.  labeled  and  scanned 

The  following  steps  determines  the  s-t  branch  connectivity  of  a 
directed  or  mixed  network: 

1.   Represent  the  network  by  a  directed  linear  graph.   Number  the 
vertices  consecutively  from  1  to  n. 
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2.  For  each  branch  (i,j)  €  p  ,  assign  c(i, j)  =  1  and  f(i,j)  =  0. 

3.  Choose  the  source  and  sink  vertices  v  and  v. . 

s      t 

4.  Apply  Algorithm  1 

Algorithm  1   (Modified  Ford  and  Fulkerson  Algorithm) 
Labeling  Routine; 

a.  Erase  all  vertex  labels 

b.  Label  v  by  (0,+).   Source  is  now  labeled  and  unscanned, 
all  other  vertices  are  unlabeled  and  unscanned. 

c.  Select  any  labeled  unscanned  vertex  v.  (including  v  ). 
Assign  labels  to  all  unlabeled  vertices  v  .  adjacent  to  v.  such  that 

1)  for  forward  branch  (i,j)  with  flow  f(i,j)  =  0,  label 
v.  by  (if+).   Vertex  v.  is  now  labeled  and  unscanned. 

2)  for  backward  branch  (j,i)  with  flow  f(j,i)  >  0,  label 

v.  by  (i,-).   Vertex  v.  is  now  labeled  and  unscanned. 
J  J 

With  all  adjacent  vertices  to  v.  labeled  and  unscanned, 

v.  is  now  scanned. 

l 

d.  Repeat  Step  c  until   either 

1)  v,    is  labeled  or 

2)  no  'more  labels  can  be  assigned  and  v,    is  unlabeled. 

In  Step  d(l),    an  augmentation  path  is  found  and  an  increase 
of  flow  is  possible.     Proceed  to  Augmentation  Routine. 

In  Step  d(2) ,    the  algorithm  terminates.      Flow  is  now  max- 
imized.    The  disconnecting  cutset  are  the  branches  joining  the  labeled 
and  unlabeled  vertices.      The  value  of  the  cut  is  the  branch  connectivity. 
Augmentation  Routine; 

a.  Let  z  =  t  and  go  to  Step  b. 

b.  If  label   on  v     is  (q,+),    set  f(q,z)   =   1.      If  label   on  v 

21  2j 

is  (qi-),    set  f(z,q)   =  0. 
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c.   If  q  =  s  return  to  Step  a  of  the  Labeling  Routine.   Otherwise 
set  z  =  q  and  go  to  Step  b  of  the  Augmentation  Routine. 
Example:   A  mixed  network 

Fig.  3.1(a)  is  the  linear  graph  representation  of  the  network. 
In  (b) ,  the  non-oriented  branch  (3 > 6)  is  split  into  two  oppositely  direct- 
ed branches  (3*6)  and  (6,3).   The  first  labeling  obtained  the  augmentation 
path  s,(s,3) >3>(3»7) »7»(?it) .  In  (c)  the  second  labeling  found  the  aug- 
mentation path  s,(s,4) ,4,(4,3) i3»(3>6) >6,(6,t) .  In  (d)  the  sink  vertex 
v,  cannot  be  labeled  since  (6,8)  and  (7«8)  are  already  saturated  and  (5»8) 


is  a  backward  branch  with  zero  flow.   The  algorithm  terminates.   The  cut 
is  |  (6,8) ,(7,8)  r    and  s-t  connectivity  is  two. 

Flow  chart  for  this  analysis  is  in  Appendix  A-l  and  the  computer 
program  MAXC0N-1A  is  in  Appendix  B-2.   To  minimize  computer  storage,  a 
vertex  label  is  further  modi.fi ed  by  using  a  single  signed  integer  (+i) 
instead  of  the  ordered  pair  (i,+).   When  i  is  negative  the  magnitude  is 
obtained  by  just  multiplying  by  minus  1 . 

In  the  program,  vertex  scanning  is  accomplished  by  using  two  stor- 
age vectors.  The  first  vector  is  a  list  of  the  labeled  vertices  being 
scanned  while  the  second  provides  storage  for  the  vertices  being  labeled. 
Once  the  scanning  process  of  all  the  vertices  listed  in  the  first  vector 
is  completed,  the  contents  of  the  second  is  transferred  to  the  first  which 
in  turn  becomes  the  list  of  vertices  to  be  scanned. 

Algorithm  1  is  applied  to  the  analysis  of  undirected  network  with 
slight  modification.  For  a  non-oriented  net,  each  branch  flow  is  con- 
sidered bidirectional  with  a  maximum  capacity  of  unity.   However,  in  no 
case  can  a  positive  back-flow  occur  since  each  augmentation  path  has  as 

its  elements  all  forward  directed  branches  with  respect  to  v  and  v. . 

s      t 
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Consequently  there  cannot  be  any  negatively  labeled  vertex.  Therefore  the 
tests  for  a  backward  branch  and  a  negative-labeled  vertex,  and  the  adjust- 
ment of  backward  flow  are  not  necessary  in  this  case. 

These  steps  are  deleted  in  Algorithm  1   and  MAXC0N-2A  (Appendix 
B-^-)   is  developed  for  determining  s-t  branch  connectivity  of  an  undirected 
network. 
3.2.2     Multi terminal  Branch  Connectivity 

For  an  n-station  network,   its  graph  would  have  n(n-l)/2  possible 
vertex  pairs.      To  apply  Algorithm  1  n(n-l)/2  times  in  order  to  obtain  the 
connectivity  between  all  the  possible  vertex-pairs  would  be  lengthy. 

Gomory  and  Hu  [_13~]  have  shown  that  for  an  undirected  network  with 
n  vertices,   any  flow  in  the  network  is  numerically  equal  to  some  flow  in 
a  maximal   spanning  tree.      Since  there  are  only  n-1   branches  in  a  spanning 
tree,    then  there  can  only  be  n-1  numerically  different  flows  possible. 
Therefore  all   the  n(n-l)/2  possible  maximal  flows  that  can  exist  in  the 
network  can  be  deduced  after  doing  n-1  flow  maximizations. 

Applying  the  above  principle  and  the  max-flow  min-cut  theorem  in 
the  analysis  of  non-oriented  nets,    the  multi-terminal  network  connectivity 
problem  can  be  solved  by  only  n-1   separate  applications  of  Algorithm  1. 
The  analysis  procedure  is  embodied  in  the  following  algorithm. 

Algorithm  2 

1.  Choose  two  vertices  v  and  v,  in  the  undirected  graph  G. 
With  branch  capacities  set  to  unity,  find  an  initial  s-t  branch  cut  using 
Algorithm  1.  Represent  this  cut  by  a  generalized  vertex  tree  with  con- 
densed vertices  X  and  X  connected  by  a  branch  with  capacity  equal  to  the 
computed  cut  value  T   . . 

2.  Ifcraw  a  condensed  network  graph  G.  where  X  is  represented 
by  a  single  vertex.   Select  two  vertices  v,  and  v,  in  X  and  find  the 
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smallest  k-1  cut  "by  Algorithm  1.   Form  the  generalized  vertex  tree  by 
dividing  X  into  X  and  X  with  *f         as  the  capacity  of  the  branch  con- 
necting them.   In  this  tree,  X  is  adjacent  to  X.  if  it  is  in  the  same 
side  of  the  cut  as  X.  ,  or  X  is  adjacent  to  X  if  it  is  in  the  same  side 

of  the  cut  as  X. . 

3.   The  above  procedure  is  repeated.   With  each  step,  a  tree 
of  generalized  vertices  connected  by  branches  with  corresponding  cut 
values  is  obtained.  To  proceed  with  the  computation, 

a.  Select  a  generalized  vertex  X. 

1 

b.  Condense  each  component  in  G  except  X.  into  a  single 
vertex  forming  the  condensed  network  graph  G.  . 

c.  Choose  two  original  vertices  v  and  v,  in  X.  as  source 

S  v  1 

and  sink  and  apply  Algorithm  1  to  G.  to  get  f     ,  . 

The  cut  obtained  in  this  process  splits  X.  into  two  parts 

X.   and  X.  ,   This  is  represented  in  the  generalized  vertex  tree  by  re- 

placing  X.  by  the  condensed  vertices  X.   and  X.   connected  by  a  branch 

1  1      X2 

bearing  the  cut  value.   All  other  branches  and  vertices  in  the  tree  are 

unchanged  except  those  which  were  formerly  connected  to  X.  .   Apply  the 
adjacency  rule  as  in  Step  2. 

k.     Repeat  the  process  until  these  generalized  vertices  con- 
sist of  exactly  one  vertex  each.   This  results  in  a  spanning  tree.   Thus 
for  any  chosen  v  and  v, ,  the  s-t  connectivity  is  the  smallest  branch 

S  L> 

weight  in  the  unique  s-t  path  in  the  tree. 

An  example  for  the  application  of  the  above  procedure  is  given 
in  Fig,  3*2.   The  generalized  vertex  tree  and  the  condensed  network  G. 
resulting  in  each  flow  calculation  are  indicated.  For  the  12-vertex  net- 
work, the  spanning  tree  is  formed  after  11  max-flow  computations. 
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A  flowchart  for  this  analysis  procedure  is  given  in  Appendix  A-2. 
This  is  further  encoded  into  a  computer  program  entitled  MAXCON-2  appear- 
ing in  Appendix  B-3.      In  this  program  the  printout  is  not  a  spanning  tree 
but  rather  a  symmetrical  branch-connectivity  matrix  which  is  essentially 
a  terminal  capacity  matrix.      Its  entries  are  derived  from  the  cut  values 
appearing  as  branch  weights  of  the  trees  generated  after  every  max-flow 
calculation  where  a  cut  isolates  a  single  vertex  from  the  rest  of  the  con- 
densed network  graph  G. .      The  connectivity  between  this  single  generalized 
vertex  and  all   the  rest  of  the  vertices  in  G  are  the  updating  entries. 

Other  information  relevant  to  MAXCON-2  are  stated  in  the  program 
documentation. 

The  concept  of  spanning  tree  cannot  be  applied  to  the  analysis  of 
a  directed  or  mixed  network  as  generally  these  networks  are  not  symmet- 
trical.     To  obtain  a  multiterminal  branch-connectivity  matrix  for  an  ori- 
ented net,   it  is  necessary  to  do  n(n-l)/2  flow  maximizations  by  Algorithm 
1.      Up  to  present,   there  is  no  available  method  that  could  reduce  this 
number  of  calculations  however  Frisch  has  developed  the  Flow  Variation 
Algorithm  [l2^   which  effectively  reduces  the  number  of  steps  required  to 

maximize  f      ,  . 
s,t 

•The  ideas  underlying  this  new  algorithm  are  as  follows:      let  v   , 

S 

v  ,  ,  and  v,  be  three  arbitrary  vertices  in  a  graph  G.  Rather  than  apply 

S  Xf 

Algorithm  1  to  independently  maximize  f  l+  and  f   , ,  first  maximize  f  ,  , 

by  Algorithm  1  then  maximize  f   ,  in  the  usual  manner  until  a  vertex  v 

s ,  t>  g 

is  labeled  for  which  f(g,V)  >  0.   Then  use  the  flow  pattern  D(s' ,t)  to 
find  F(s,t). 

Applying  this  basic  idea  to  solve  the  multiterminal  connectivity 
problem  which  is  a  special  case  of  a  general  flow  problem,  all  that  is 
necessary  is  to  slightly  modify  the  algorithm  to  fit  into  the  problem. 
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The  result  is  a  simplified  version  of  the  original  algorithm.     The  following 

notations  are  used: 

TT      .        -  flow  path  for  vertex-pair  v     and  v, 
s,t  *  *  s  t 

TT   ,    ,      -  flow  path  for  vertex-pair  v   ,    and  v. 

F(s,t)    -  flow  pattern  for  vertex-pair  v     and  v, 

s  x» 

D(s',t)-  flow  pattern  for  vertex-pair  v   ,    and  v, 

s  x> 

f(i,j)   -  a  branch  flow  in  path  n     . 

d(i,j)   -  a  branch  flow  in  path  n   ,    , 

S    i  x» 

tt  -  section  of  path  rr     ,    from  v     to  v     called  "Augmenta- 

s,g  s,t  s  g  & 

tion  Path" 

T 
n  ,    -  section  of  path  n  ,  ,  from  v  ,  to  v  called  "Trunca- 
s' ,g  *  s' ,t       s'     g 

tion  Path" 

u 

tt   .    -  section  of  path  n  ,  ,  from  v  to  v,  called  "Homing 

g»  U  S  | "t  g       t/ 

Path" 

u 

L      -  a  loop  detected  by  the  Homing  Routine;  m  is  the  first 

m,m         r  *  ° 

vertex  in  the  loop 

T 
L      -  a  loop  detected  by  the  Truncation  Routine;  p  is  the 
P.P      . 

first  vertex  in  the  loop 
The  subscript  c  whenever  it  appears  with  the  symbols  indicates  the 
c   stage  of  max-flow  calculation. 

The  following  routines  are  employed  in  the  algorithm: 
Routine  A  -  Detecting  an  Augmentation  Path 

1.  Label  v  by  (s,+).   v  is  now  labeled  and  unscanned.   All 

s  s 

other  vertices  are  unlabeled  and  unscanned. 

2.  Select  any  labeled  and  unscanned  vertex  v.  .   Assign  labels 
to  all  unlabeled  vertices  adjacent  to  v.  such  that 
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a.  for  forward  branch  (i,j)  €  P   with  flow  f(i,j)  =  0, 
label  v.  by  (i,+),   v.  is  now  labeled  and  unscanned. 

b.  for  backward  branch  (j,i)  6  P   with  flow  f(j,i)  >  0, 
label  v.  by  (i,-).   v.  is  now  labeled  and  unscanned. 

With  all  adjacent  vertices  to  v.  labeled  and  unscanned,  v. 
is  now  scanned, 

3.  Repeat  Step  2  until  either: 

a.  v,  is  labeled  or 

b.  a  vertex  v     is  labeled  with  d  (g,V)   >  0 

g  cVb 

c.  no  more  vertices  can  be  labeled  and  v,  is  unlabeled. 
Observe  that  except  for  Step  3(b)  this  routine  is  the  Ford 

and  Fulkerson  Labeling  Routine. 

Routine  T  -  Detecting  a  Truncation  Path 
1.  Let  z  =  g 


2.  Find  one  vertex  v.  such  that  d  (i,z)  >  0  and  label  v.  by 

1  c  1 


(■) 


3.   Let  z  =   i 

a.  If  v  has  assumed  a  value  i  before,  terminate. 

z 

b.  If  i  =  s' ,  terminate. 

c.  Otherwise  return  to  Step  2. 

T 
If  the  routine  terminates  in  Step  3(a),  a  loop  L    has  been 

p,p 

detected  rather  than  a  path. 

Routine  H  -  Detecting  a  Homing  Path 
1.  Let  q  =  g 


(q) 


2.   Find  one  vertex  v.  such  that  d(q,j)  >  0  and  label  v.  by 


3.   Let  q  =  j 

a.   If  q  has  assumed  the  value  j  before,  terminate. 
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b,  If  j  =  t,  terminate. 

c.  Others se  return  to  Step  2. 

H 
If  the  routine  terminates  in  Step  3(a),  a  loop  L    has  been 

'         m,m 

detected  instead  of  a  path. 

The  main  algorithm  for  obtaining  the  multiterminal  branch- connec- 
tivity of  an  oriented  network  is  now  stated  in  terms  of  Routines  A,  T, 
and  H. 

Algorithm  3 

1.  Choose  two  vertices  v  and  v..   Using  Routine  A,  find  an 

s  x> 

initial   s-t  cut. 

A 

2.  Using  Routine  A  find  an  s-g  augmentation  path  tt         or  an 

s  >  g 

A 
s-t  augmentation  path  n   .  In  the  former  case  go  to  Step  3.  In  the  lat- 

s  $  x> 


ter  case  go  to  Step  8.   If  no  such  path  can  be  found,  terminate. 

H 

T 

e 


3.   Using  Routine  H  find  a  Homing  Path  tt   ,  or  a  Homing  Loop 

g>  t 


1_r 

L   .  In  the  former  case  go  to  Step  k;   in  the  latter  case  go  to  Step 
m,m 

5(a). 

T 
4.      Using  Routine  T  find  a  Truncation  Path  tt   ,    ,    or  a  Trunca- 

S   J  u 

T 
tion  Loop  L   .   In  the  former  case  go  to  Step  6;  in  the  latter  case  go 

P*P 
to  Step  5(b). 

5«   a.   Set  to  zero  all  d  (i,j)  for  all  branches  (i,j)  in  the 

loop  L   .   Erase  all  vertex  labels  and  go  to  2. 


m,m 


b.   Set  to  zero  all  d  (i,j)  for  all  branches  (i,j)  in  the 


T 
loop  L   .   Erase  all  vertex  labels  and  go  to  2. 


P»P 


H 


6.   Trace  and  erase  g-t  path.   For  every  branch  (i,j)  in  TT   ,  , 


set 


fc+1(i.o)  -  1 
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7.  Truncate  the  s'-g  path.   For  every  branch  (i,j)  in  n  ,   , 

s  ,  g 

set  dc+1(i,j)  =  0. 

8.  Augment  the  s-g  path.  For  every  branch  (i,j)  in  tt   ,  set 

s»g 

fc+1(i.a)  -  i. 

Erase  vertex  labels  and  return  to  Step  2. 

If  a  "loop  of  flow"  is  detected,  the  flow  in  this  loop  must  be 

set  to  zero  as  in  Step  5  as  otherwise  a  nonfeasible  flow  will  be  achieved. 

As  an  example,  consider  the  graph  in  Fig.  3«3.   Let  the  indicated 

flow  pattern  D(s' ,t)  be  the  result  of  initial  maximization  of  f  ,  , . 

S  f  x» 

Designate  v„  as  the  source  for  the  second  stage  of  maximization.   Step  2 

of  the  algorithm  labels  v_  and  detects  it  as  a  vertex  v  .   Step  3  finds 

3  g 

the  sequence  3>  (3>6)  ,6,(6,8)  ,8  as  the  Homing  Path  tt  ,.  The  path  s',(l,3)i 

Si  "& 

T  A 

3  is  the  Truncation  Path  tt   ,      ,    and  s,(2,3)>3  the  Augmentation  Path  tt 

s  , g  s,g 

Step  6  sets  to  zero  the  branch  flows  d  (i,j)  from  v~  to  v,  and  replaces 

C  5       T/ 

each  branch  flow  in  the  path  by  f^.1(i»j)  =  it   Step  7  truncates  s'-3  path 

by  setting  d  (1,3)  to  zero,  and  Step  8  augments  s-3  path  by  setting 

f  .(2,3)   !■   One  s-t  path  is  now  processed  and  the  algorithm  will 

search  foj   uother  one  until  no  more  augmentation  paths  can  be  found.  By 

then,  f   ,  will  be  maximal, 
s,  t 

.  Flowchart  of  this  algorithm  is  in  Appendix  A- 3  and  its  computer 
implementation  MAXCON-1  is  presented  in  Appendix  B-l. 


Fig.  3-3 
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3.3  ANALYSIS  WITH  RESPECT  TO  VERTEX  DISCONNECTION 

The  criterion  of  survivability  attached  to  this  problem  is  the  number 
of  stations  that  must  be  deleted  from  the  network  in  order  to  disrupt  com- 
munication.  This  is  the  smallest  disconnecting  set  W,  of  a  network 
graph.   Given  any  pair  of  vertices  v  and  v, ,  assign  unit  capacities  to 
all  vertices  except  v  and  v,  (which  have  infinite  capacities),  then  the 
smallest  s-t  vertex  cutset  Oi  is  numerically  equal  to  the  maximum  num- 

ber  of  vertex  disjoint  s-t  paths.   The  value  of  (J    is  the  network  ver- 

S  ,  *C 

tex-connectivity  with  respect  to  the  given  source  and  sink. 
3.3-1   s-t  Vertex  Connectivity 

To  solve  the  problem  of  finding  W   , ,  each  vertex  in  G  is  split 
into  two  vertices  connected  by  a  branch  of  unit  capacity  forming  a  new 
graph  G'.   Algorithm  1  is  then  applied  to  G'.   The  result  is  a  branch 
cutset  composed  only  of  the  unit  capacity  branches.  The  regular  network 
branches  are  considered  having  infinite  capacities.   This  approach  is 
however  very  inefficient  since  the  number  of  vertices  are  almost  doubled. 

An  algorithm  known  as  Vertex-Pair  Connectivity  derived  from  Algo- 
rithm 1  was  developed  by  Frisch  [ll^j  which  applies  directly  to  graphs  of 
directed,  undirected,  or  mixed  networks  without  splitting  the  vertices. 
The  vertex-splitting  process  is  absorbed  in  the  modification  of  the  label- 
ing routine  of  Algorithm  1.   The  algorithm  is  stated  here  in  a  different 
form.   Several  terms  and  rules  must  first  be  defined: 

1.  A  vertex  is  labeled  with  an  ordered  triple  (i,+,M).   The 
first  and  second  entries  are  defined  as  in  the  original  Labeling  Routine. 
The  third  entry  is  a  zero  or  one  as  will  be  explained  in  the  algorithm. 

2.  A  vertex  is  strongly  labeled  if  the  second  entry  is  a  minus, 
and  weakly  labeled  if  it  is  a  plus. 

3.  A  circle  around  the  second  entry  denotes  the  vertex  is 

scanned. 
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Ur.      If  a  vertex  is  weakly  labeled  and  unscanned,  and  then 
receives  a  strong  label,  the  weak  label  is  erased.   However  if  the  label 
is  scanned,  it  should  not  be  erased. 

Algorithm  h 

Labeling  Routine 

1.  Erase  vertex  labels.   Label  v  by  (s,+,0).   v  is  now 

s  s 

labeled  and  unscanned.   All  other  vertices  are  unlabeled  and  unscanned. 

2.  Select  any  labeled  unscanned  vertex  v.  .   If  none,  ter- 
minate. 

a.   Suppose  v.  is  weakly  labeled: 

1)  If  f(i,j)  =  0,  (i,j)  6  P  and  j  ^  s,  weakly 
label  all  vertex  v.  by  (i,+,0).   v.  is  now  weakly  labeled  and  unscanned. 

2)  If  f(j,i)  -  1,  (j,i)  6  r  ,  j  /  s,  and  v  not 
strongly  labeled,  strongly  label  v.  by  (i,-,0).  v.  is  now  strongly  label- 
ed  and  unscanned. 

Scan  v.  by  encircling  the  +  in  the  v.  label, 
b.   Suppose  v.  is  strongly  labeled 

1)  If  f(i,j)  =  0,  (i,j)  €  T  ,  and  j  /  s,  weakly 
label  all  unlabeled  vertices  v.  by  (i,+,M).  v.  is  now  weakly  labeled  and 
unscanned. 

2)  If  f(j,i)  =  1,  (j,i)  €  V  ,  j  ^  s,  and  v  not 
strongly  labeled,  strongly  label  v.  by  (i,-,M).  If  v.  is  weakly  labeled, 
erase  the  weak  label,  v.  is  now  strongly  labeled  and  unscanned. 

In  (b),  if  v.  is  both  weakly  and  strongly  labeled, 
set  M  =  lf  otherwise  M  =  0.  Scan  v.  by  encircling  the  minus  sign  in  the 
v.  label, 

3.  If  v,  is  labeled,  proceed  to  Augmentation  Routine. 


Otherwise  return  to  Step  2. 
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Augmentation  Routine: 

1.  Let   z  =  t 

2.  a.   If  v  is  weakly  labeled  by  (q,+,M),  set  f(q,z)  =  1. 
b.   If  v  is  strongly  labeled  by  (q,-,M),  set  f(z,q)  =  0. 

3.  If  q  =  s  return  to  Step  1  of  Labeling  Routine.   Otherwise 
proceed  to  Step  k. 

4.  Let  eT=   q.  If  M  =  1  return  to  Step  2(b).   Otherwise  return 
to  Step  2(a) o 

To  identify  the  s-t  vertex  cutset,  let  U  be  the  set  of  all  unlabeled 

vertices,  and  L  be  the  set  of  all  labeled  vertices.   Further  let  L  be  the 

u 

set  of  labeled  vertices  connected  to  an  unlabeled  vertex  by  a  branch  di- 
rected from  a  labeled  vertex  to  an  unlabeled  one,  then  the  s-t  vertex  cut 

X  is  given  by: 

m    &      J 

Xm  =  Lu  if  s  /  Lu  (3.3.1.1) 

Xm  =  (Lu  "  VS)  U  {vh:(s'h)  €   (L'U))  lf  s  €  Lu  (3.3.1.2) 

Since  the  source  vertex  cannot  be  a  component  of  the  vertex  cut,  Eq. 

(3.3.1.2)  shows  that  if  v  is  connected  to  unlabeled  vertices  v,  ,  then 
N       '  s  h 

the  vertices  v,  are  some  component  of  the  cut. 

The  graph  in  Fig.  3«^  illustrates  the  application  of  the  algorithm 
for  a  directed  network.   Assume  that  the  path  indicated  by  solid  arrows 
is  a  result  of  initial  labeling  and  augmentation.   Assume  further  that 
the  second  application  of  labeling  routine  is  just  completed  and  v.  is 
labeled  as  indicated  by  the  vertex  labels.   When  the  augmentation  routine 
is  applied,  it  will  backtrack  the  path  in  the  usual  manner,  however,  on 
reaching  v^  it  has  to  select  the  strong  label  since  v,   has  M  =  1.   If  the 
weak  label  were  selected,  a  non-feasible  flow  would  result.   This  explains 
the  purpose  of  the  third  entry  in  the  label.   The  only  case  where  two 


3^ 


(4,0,0) 
f3 


(sAq 


8)^,0,0) 


s,+,o) 


Fig.    3.4 
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scanned  labels  can  occur  on  a  vertex  is  the  situation  in  v.  where  a  flow 
of  unit  magnitude  already  exists  through  the  vertex  and  this  vertex  re- 
ceives labels  in  the  succeeding  iteration  of  the  algorithm.   If  the  weak 
label  in  v^  was  unscanned  when  it  received  a  strong  label  then  the  weak 
label  should  have  been  erased  and  v^  would  have  been  labeled  (^,+,0). 

Program  MAXC0N-3A  in  Appendix  B-6  is  a  computer  implementation  of 
Algorithm  4,   A  different  way  of  circumventing  the  situation  in  vertex 
Vj,  of  Fig.  J>,k   was  applied.  For  a  vertex  like  v^  where  a  flow  already 
exists,  if  in  the  succeeding  labeling  process  v.  receives  a  weak  label, 
it  will  not  be  scanned,  but  if  it  receives  a  strong  label  it  will  be  scan- 
ned.  The  latter  would  indicate  that  the  existing  outward  flow  from  v^   can 
in  effect  be  diverted  to  another  path.  In  the  former,  it  would  be  waste- 
ful to  continue  seeking  another  path  passing  through  v^  when  the  vertex 
is  already  saturated.   In  no  case  then  can  a  vertex  acquire  two  scanned 
labels.  The  third  entry  in  the  label  is  therefore  no  longer  necessary. 

The  same  form  of  label  as  in  MAXC0N-1A  is  used  in  this  program. 
All  steps  in  the  algorithm  that  are  associated  with  the  third  entry  in 
the  label  are  deleted. 

A  flowchart  of  .Algorithm  4  is  in  Appendix  A-4. 
3.3,2  Multiterminal  Vertex  Connectivity 

The  Flow  Variation  Algorithm  of  Frisch  as  presented  in  Sec.  3»2 
maximizes  flows  in  a  finitely  branch- weighted  network,  and  it  was  readily 
adapted  to  solving  the  multiterminal  branch-connectivity.  In  this  sec- 
tion the  problem  is  different  in  the  sense  that  the  vertices  control  the 
information  flow. 

In  order  to  utilize  the  principle  of  the  Flow  Variation  Algorithm 
here,  the  Routine  A  of  the  algorithm  must  be  altered.   Recall  that  Routine 
A  is  the  Ford-Fulkerson  Labeling  Routine.   If  this  routine  is  replaced  by 
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the  Labeling  Routine  of  Algorithm  k  which  is  Frisch's  Vertex-Pair 
Connectivity  Algorithm,   it  would  enable  the  Flow  Variation  Algorithm  to 
handle  the  multiterminal  vertex-connectivity  problem.      Consequently  the 
method  in  which  the  s-g  and  s-t  paths  are  augmented  must  also  conform 
with  the  Augmentation  Routine  of  Algorithm  4. 

With  these  changes,   the  algorithm  flowchart  is  the  same  as  in  Ap- 
pendix A-3,    however  there  is  a  significant  change  in  the  computer  pro- 
gram MAXCON-3  as  presented  in  Appendix  B-5  compared  to  MAXCON-i. 

The  solution  given  earlier  for  obtaining  the  multiterminal  vertex- 
connectivity   of  a  directed  network  applies   equally  well  for  undirected 
networks.      Sometimes  one  finds  it  impractical   to  compute  all  the  entries 
in  the  entire  connectivity  matrix  especially  when  investigating  a  very 
large  network.      It  is  often  sufficient  to  just  verify  whether  a  network 
satisfies  a  certain  minimum  degree  of  connectivity. 

With  this  analysis  criterion  it  is  only  necessary  to  verify  that 
the  network  is  at  least  R-connected.  The  connectivity  being  referred  to 
here  is  the  vert ex- connectivity  which  we  know  is  the  smallest. 

The  following  theorems  of  Kleitman  [l^]  greatly  reduce  the  number 
of  verifications  necessary  rather  than  doing  n(n-l)/2  separate  tests  for 
the  existence  of  R  vertex-disjoint  paths. 

Theorem: 

To  verify  that  all  n(n-l)/2  vertex  pairs  in  G  can  be  connected 
by  at  least  R  vert ex- disjoint  paths,  it  is  only  necessary  to  verify  that 
some  set  of  R  vertices  are  each  connected  by  at  least  R  vert ex- disjoint 
paths  with  every  vertex  in  the  graph. 

Two  other  theorems  resulting  from  the  above  are  stated  here  in  an 
algorithm  form  and  are  in  themselves  the  analysis  procedure. 
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Algorithm  5 


1.  Choose  any  vertex,  say  v.  as  source. 
Let  i  =  1. 

a.      Select  a  vertex  v.  in  G.    as  sink.      Remove  from  G.    all 

J     i  i 

r        \ 


'     ,   k  =  1,...,1  previously  designated  as  sink  that 


other  vertices  •  v. 

are  adjacent  to  v..   If  1  >  R  go  to  Step  l(b). 

•J 

J.)  For  v.  connected  to  v. 
J  i 

Excluding  the  direct  path  i-j,  use  Algorithm  4  to 
verify  existence  of  q  vert  ex- disjoint  paths  from  v.  to  v..   Set  s  =  q+1. 
If  s  <  R  go  to  Step  4.  Otherwise  go  to  Step  1(b). 

b.  Are  all  n-i  vertices  processed?  If  not  go  to  Step  1(a), 
otherwise  go  to  Step  2. 

2.  Decrease  R  by  i.   If  R  =  0,  terminate;  otherwise  go  to 
Step  3. 

3.  Remove  v.  from  G. .  Choose  vertex  v.  .  as  source.   Increase 

11  l+l 

i  by  one.      Return  to  Step  1(a). 

4.  Record  s-connectivity  between  v.    and  v..      Terminate.      Net- 
work is  not  R-connected. 

The  first  theorem  may  be  restated  in  a  more  specific  sense: 

If  there  are  at  least  R  vert ex- disjoint  paths  between  v.  and 
v.  and  v.  and  v.  in  G,  then  to  verify  the  existence  of  R  such  paths  be- 
tween v.    and  v.  in  G,   we  need  only  look  for  R  -  1   such  paths  in  G  -  iv.l    . 

For  the  proof  of  these  theorems  and  illustrative  example,    the  read- 
er is  referred  to  the  original  paper  by  Kleitman. 

A  flowchart  for  computer  implementation  is  given  in  Appendix  A-5. 
Program  MAXCON-4  in  Appendix  B-7  is  encoded  to  analyze  an  undirected  net- 
work for  a  uniform  R- connectivity  where  the  value  of  R  is  initially 
specified. 
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IV.      DESIGN  OF  LOW-COST  SURVIVABLE  NETWORK 

4.1  INTRODUCTION 

There  are  now  a  number  of  synthesis  procedures  for  survivable  networks 
appearing  in  literatures  elsewhere.  These  are  summarized  by  Frank  and 
Frisch  in  [l0[]  and  are  categorized  according  to  link  or  station  damage. 
The  synthesis  approach  are  all  aimed  towards  maximization  of  redundancy 
between  stations  and  satisfying  network  redundancy  requirement  with  the 
assumption  that  cost  is  either  uniform  for  each  link/station  or  not  con- 
sidered at  all. 

In  the  design  of  practical  flow  networks,  the  cost  of  construction 
increases  with  increase  in  redundancy  therefore  the  design  approach  must 
take  into  account  non-uniform  cost  for  the  construction  of  links  or  sta- 
tions.  The  goal  is  a  minimum  cost  network  that  satisfies  a  prespecified 
redundancy  requirement.  So  far  there  is  no  exact  analytical  method  that 
could  be  applied  to  solve  the  problem.   A  procedure  to  search  for  a  low- 
cost  network  with  respect  to  vertex  disconnection  is  presented  by  Steig- 
litz,  etal  [_1Q~]   using  an  heuristic  approach.  The  algorithms  they  devel- 
oped were  based  on  the  techniques  for  an  approximate  solution  to  the 
traveling  salesman  problem.   The  salesman's  tour  is  a  Hamiltonian  circuit 
which  is  a  special  case  of  the  flow  network  having  a  redundancy  of  two. 

4.2  DESIGN  WITH  BRANCH  DISCONNECTION  CRITERION 

In  some  cases  branch  destruction  is  more  likely  than  vertex  destruc- 
tion and  it  is  just  as  necessary  to  obtain  a  design  method  using  the 
former  criterion  as  the  latter.   A  procedure  is  presented  in  this  chapter 
for  design  of  a  low  cost  network  (with  high  probability  of  being  optimal) 
with  at  least  r.  .  branch  disjoint  paths. 
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The  same  procedure  used  by  Steiglitz,  etal  is  employed  consisting  of 
two  parts.   The  first  part  called  the  Starting  Routine  generates  a  fea- 
sible initial  network  (one  which  satisfies  redundancy  requirement)  and 
the  second  part  is  the  Optimizing  Routine  which  makes  an  iterative  im- 
provement on  the  network  cost  "by  branch  exchanges  while  maintaining  fea- 
sibility. Portions  of  the  procedure  that  differ  from  that  of  Steiglitz 
are  the  feasibility  testing  of  the  network  and  the  optimizing  routine. 
The  method  of  optimization  that  will  be  presented  is  an  adaptation  of  the 
recently  published  heuristic  algorithm  for  the  traveling  salesman  problem 
by  Lin  and  Kernighan  [l6]]. 

The  heuristic  design  algorithm  is  stated  as  follows: 

Algorithm  6 

1.  Generate  a  feasible  initial  network  G  that  satisfies  the 
specified  redundancy  R  using  the  Starting  Routine. 

2.  Attempt  to  find  an  improved  feasible  solution  G*  by  some 
local  transformation  using  the  Optimizing  Routine, 

3.  If  a  similar  network  with  a  lower  cost  is  found,  then  re- 
place G  by  G'  and  repeat  from  Step  2. 

k.     If  no  local  improvement  can  be  found,  G  is  a  locally  optimal 
solution.  Repeat  from  Step  1  until  computation  time  runs  out  or  the  so- 
lution is  satisfactory. 
Jj-,2.1  The  Starting  Routine 

This  routine  was  developed  by  Steiglitz,  etal  in  [.18]]. 

1.  Plot  vertices  representing  stations  in  a  Euclidian  2-space. 

2.  Number  the  vertices  uniformly  at  random. 

3.  Assume  cost  to  be  proportional  to  distance  and  obtain  the 
cost  matrix  from  the  Euclidian  distances  between  vertices  truncated  to  an 
integer. 


4.  Assign  redundancy  to  each  vertex. 

5.  Add  branches  to  the  network  one  at  a  time.   Reduce  "by  unity 
the  assigned  vertex  requirement  every  time  a  link  is  connected  to  that 
particular  vertex,  A  branch  is  added  in  this  manner: 

a.  Connect  the  vertex  with  highest  updated  requirement  to 
another  vertex  with  the  next  highest  updated  requirement;  for  all  vertices 
having  equal  requirement,  choose  one  that  results  in  lower  cost. 

b.  In  cases  of  ties  in  requirement  and  cost,  these  are 
resolved  by  choosing  the  vertex  highest  on  the  list. 

c.  No  parallel  branches  are  allowed. 

6.  When  the  requirement  in  all  the  vertices  are  reduced  to 
zero,  test  feasibility  of  the  initial  solution  using  Algorithm  2. 

7.  If  solution  is  feasible  go  to  Optimizing  Routine,  otherwise 
repeat  from  Step  2. 

This  routine  is  made  nondeterministic  by  numbering  the  vertices 
uniformly  at  random  at  each  start.  Randomized  starting  solution  has  been 
observed  experimentally  to  produce  a  variety  of  starting  networks  as  con- 
trasted to  constructive  solutions  which  are  deterministic  where  more  than 
one  solution  may  not  be  obtained. 

The  analysis  tools  presented  in  Chapter  III  come  in  handy  for  test- 
ing feasibility.  For  the  initial  network,  the  one  that  fits  is  Algorithm 
2  which  determines  the  overall  branch  connectivity  of  the  network.   A 
complete  checkout  need  not  be  carried  any  further  once  inf easibility  is 
detected. 
4.2.2  Optimizing  Routine 

The  optimizing  procedure  applied  by  Steiglitz  is  the  X-change  where 
two  branches  are  interchanged  at  a  time,  and  determined  whether  the  trans- 
formation is  favorable  and  preserves  feasibility.  Lin  and  Kemighan  have 
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lately  experimented  on  a  more  systematic  and  effective  way  of  branch 
transformation  in  a  network  with  2-redundancy.   Their  results  are  applied 
here  for  a  network  with  redundancy  greater  than  two. 

The  basic  ideas  underlying  the  transformation  are  as  follows: 
Consider  a  nonoptimal  but  feasible  network  G  with  cost  c(G)  and  an  opti- 
mal or  near  optimal  network  G'  with  the  same  number  of  vertices  and  with 
cost  c(G*)  <  c(G) .  G  is  nonoptimal  because  it  has  k  links  x  ,  ...,x,  that 
are  "out  of  place".  To  make  G  optimal  they  'should  be  replaced  by  k  links 
y1i...»yk  of  G' .   We  try  to  find  k  by  identifying  sequentially,  element 
by  element,  two  disjoint  sets  of  links  X  =  \   x,  , . . .  ,  x,  I  and  Y  =  <  y .  , . . ,  ,y  [ 
using  a  gain  criterion.   Let  the  costs  of  x.  and  y.  be  I  x.  I   and  ly.l 
respectively.   Define  g.  =  I  x.  I  -  I  y.  I   as  the  gain  from  exchanging  x. 

with  y.  .   For  some  i,  g.  may  be  negative  however  we  want  the  final  gain 

1  k1 

to  be  positive,  i.e.,  Eg.  >  0.   We  are  looking  for  sequences  of  g. ' s  whose 

i  X  k  1 

partial  sum  is  always  positive.   When  Eg.  <  0,  the  search  is  stopped. 

i 
The  k  links  in  X  are  deleted  or  "broken"  and  replaced  by  k  links 

in  Y  resulting  in  G*  which  has  lower  cost.  The  process  is  iterated  until 

no  further  reduction  in  cost  is  attained. 

We  now  state  the  routine  more  formally  assuming  that  a  starting 

solution  has  been  generated: 

1.  Set  P*  =  0.   (P*  is  the  best  improvement  so  far).   Choose 
any  vertex  as  t.  and  let  x,  be  one  of  the  edges  of  G  incident  to  t. .   Let 
i  =  1. 

2.  Designate  the  other  endpoint  of  x.  as  t  .   From  t  choose 
y.  different  from  existing  links  in  G  such  that  g.X)  and  designate  the 
other  endpoint  of  y.  as  t~.   If  no  such  y,  exists,  go  to  Step  5(d). 

3.  Let  i  =  i  +  1,  Choose  x.  (which  currently  joins  t   .  to 
t_.  )  and  y.  as  follows: 

42 


a.  x.  is  chosen  so  that  it  is  the  longest  link  in  G 

incident  to  vertex  t_.  , 

2i-l 

b.  y.  is  some  link  that  can  be  drawn  from  vertex  t  .  sub- 
ject to  (c),  (d),  and  (e)  below.   If  no  y.  exists,  go  to  Step  k, 

c.  To  guarantee  that  >:'  s  and  y's  are  disjoint,  x.  cannot 

be  a  link  previously  joined  (i.e.,  ay.,  j  <  i) ,  and  similarly  y.  cannot 

be  a  link  previously  broken. 

i 


d.  P.    =  £  g.  >  0.      (Gain  criterion) 

i   1  J 

e.  y.  chosen  must  permit  breaking  of  an  x.  , 


f.   Before  y.  is  drawn,  check  if  closing  up  by  joining 

t  .  to  t.  will  give  a  gain  value  better  than  the  best  obtained  previously, 

Let  y  *  be  a  link  connecting  t?.  to  t,  and  let  g.  *  =  I  x.  I  -  |y-*|   • 

If  P.  ,  +  g.*  >  P*,  set  P*  =  P.  ,  +  g.  *  and  let  k  =  i.   (P*  is  now  the 
l-l   Di  l-l   °i  v 

best  improvement  in  G.  P*  >  0  and  must  be  monotonically  nondecreasing. 
The  index  k  defines  the  sets  to  be  exchanged  to  achieve  P*) . 


either 


4.   Terminate  the  search  for  x.  and  y.  in  Steps  1  to  3  when 


a.  no  further  links  x.  and  y_.  satisfy  3(c)  to  (e). 


b.  P.  •<  P* 
i  — 

Delete  the  x. 's  and  replace  them  by  the  y. 's,  i  =  l,...,k, 
to  obtain  G' . 

5.  Test  the  feasibility  of  G*  (see  Sec.  4.2.3).  If  feasible, 
take  G'  as  the  starting  network  and  repeat  from  Step  1.  Otherwise  gen- 
erate another  initial  solution  using  the  Starting  Routine. 

6.  If  P*  =  0,  a  limited  backtracking  is  invoked: 

a.  Repeat  Steps  3  sJid.  k-,   choosing  y's  in  increasing  cost 
as  long  as  they  satisfy  the  gain  criterion  g.  +  g2  >  0. 
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b.  If  all   choices  of  y_  in  Step  3(b)    are  exhausted  without 
profit,   return  to  Step  3(a)   and  try  the  alternate  choice  for  x  „ 

c.  If  this  also  fails  to  give  improvement,    a  further  back- 
up is  performed  to  Step  2,   where  the  y/s  are  examined  in  increasing  cost. 

d.  If  the  y. ' s  are  also   exhausted  without  profit,    try  the 
alternate  x,    in  Step  1. 

e.  If  this  fails,    select  another  vertex  t,    and  repeat 
from  Step  1. 

7.      The  routine  terminates  when  all  vertices  in  G  have  been 
designated  t.    without  profit.      Return  to  the  Starting  Routine. 

A  justification  of  Step  6  is  in  order:      Experimental  results 
by  Lin  and  Kernighan  showed  that  for  the  2-redundancy  network,   backtrack- 
ing on  the  third  level  and  higher,   i.e.,   i  >  3»   indicated  considerable 
time  penalty,   hence  the  only  alternates  for  x, ,   y.    and  x   ,y     at  the  first 
and  second  levels  are  examined.      Other  refinements  to  the  procedure  dis- 
cussed in  the  original  paper  are  not  applied  here  since  actual   computer 
runs  may  indicate  different  time-cost  trade-offs  for  networks  with  redun- 
dancy greater  than  2. 
^.2.3     Feasibility  Test 

The  optimizing  routine  for  2-redundancy  network  has  a  built-in 
feasibility  test  but  is  not  applicable  to  higher  redundancy  networks. 
The  fact  that  G  is  feasible  does  not  imply  that  G1    is  also  feasible  even 
though  the  degree  of  the  vertices  satisfies  the  requirement.     It  is  there- 
fore necessary  to  test  the  redundancy  after  every  transformation.     Rather 
than  applying  Algorithm  2  to   examine  the  new  network  G' ,    the  number  of 
flow  calculations  can  be  reduced  significantly  using  the  theorems  stated 
by  Steiglitz,    etal  for  an  X-change  or  2-change.      Extending  this  to  a  k- 
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change  transformation  and  with  uniform  redundancy,  we  need  only  to  check 
the  redundancy  between  the  pairs  of  vertices  whose  connecting  links  were 
deleted.  Algorithm  1  is  most  suitable  for  investigating  feasibility  after 
every  any  transformation. 

Exampl e  s 

Given  a  set  of  stations  whose  geographical  arrangement  is  trans- 
ferred on  an  Euclidian  plane  as  in  Fig.  5«1».  w©  attempt  to  construct  a 
minimum  cost  flow  network  with  a  uniform  redundancy  of  3.  We  model  the 
network  as  a  graph  of  vertices  and  branches  and  illustrate  the  applica- 
tion of  the  design  procedure. 

Start  by  numbering  the  vertices  uniformly  at  random  as  in  Fig.  5,1. 
Measure  the  distances  between  vertices  (assuming  proportional  to  cost)  and 
obtain  the  cost  matrix.   Add  branches  one  at  a  time  and  keep  track  of  the 
number  of  incident  branches  on  each  vertex  by  using  an  updated  require- 
ment table  as  shown.   Row  1  of  this  table  is  filled  up  by  the  assigned 
vertex  requirements.   This  row  has  all  the  same  entries  so  start  with 
vertex  1.   The  cost  matrix  shows  branch  (1,6)  has  the  least  cost  (here 
the  parentheses  are  used  in  place  of  brackets  to  represent  undirected 
branches).   Encircle  this  entry  in  the  matrix  and  construct  the  branch 
(1,6)  in  Fig.  5»2.  Form  row  2  of  the  requirement  table  by  reducing  en- 
tries for  vertices  1  and  6  by  unity.  The  highest  updated  requirement  in 
this  row  is  3  and  the  highest  on  the  list  is  vertex  2.  From  the  cost  ma- 
trix we  see  that  branch  (2,7)  has  the  least  cost.   Encircle  this  entry 
and  construct  the  branch  (2,?).  The  procedure  is  continued  and  the  se- 
quence of  branches  constructed  next  are  (3»^)»  (5,8) , (l ,?) ,(2,8) ,(3,6) , 
(4,5), (1,2), (3, 7), (4,8),  and  (5,6). 
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Requirement  Table 


Vertex  List 

1 

2 

3 

4 

5 

6 

7 

8 

Branch 

Requirement 

3 

3 

3 

3 

3 

3 

3 

3 

(1.6) 

2 

3 

3 

3 

3 

2 

3 

3 

(2,7) 

2 

2 

3 

3 

3 

2 

3 

3 

(3,4) 

2 

2 

2 

2 

3 

2 

2 

3 

(5,8) 

2 

2 

2 

2 

2 

2 

2 

2 

(1.7) 

1 

2 

2 

2 

2 

2 

1 

2 

(2,8) 

1 

1 

2 

2 

2 

2 

1 

1 

(3,6) 

1 

1 

1 

2 

2 

1 

1 

1 

(4,5) 

1 

1 

1 

1 

1 

1 

1 

1 

(1.2) 

0 

0 

1 

1 

1 

1 

1 

1 

(3,7) 

0 

0 

0 

1. 

1 

1 

0 

1 

(4,8) 

0 

0 

0 

0 

1 

1 

0 

0 

(5,6) 

0 

0 

0 

0 

0 

0 

0 

0 
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When  the  number  of  vertices  is  odd,  the  last  row  of  the  requirement 
table  has  one  vertex  left  with  an  updated  requirement  of  1  while  the  rest 
are  zero.   Connect  this  vertex  to  one  which  produces  the  least  branch  cost, 
at  the  same  time  satisfying  Step  5(c)  of  the  starting  routine. 

The  initial  solution  generated  by  the  starting  routine  is  shown  by 
solid  lines  in  Fig.  5»2.  It  is  found  feasible  using  Algorithm  2.  The 
cost  of  the  network  is  obtained  by  summing  up  the  values  of  the  encircled 
entries  in  the  cost  matrix.  For  this  example,  the  initial  cost  is  261. 

We  apply  now  the  Optimizing  Routine  to  reduce  the  cost.   Set  P*  = 
0  and  i  =  1.   In  Fig.  5«2,  choose  vertex  1  as  t. .   x,  is  selected  such 
that  it  is  the  longest  among  the  three  incident  branches  to  t. .  Vertex 
2  which  is  the  other  endpoint  of  x.  is  designated  t  .   From  tr    find  a  y, 
using  Step  3  of  the  routine  such  that  g>  >  0.   Only  the  imaginary  branch 
(2,4)  is  available  since  the  other  possibilities  (2,5) » (2,6),  and  (2,3) 
would  result  to  g.  <  0.  Choosing  (2,4)  as  y.  the  gain  g.  =   J  x,  I 
jy. [   =3-   The  overall  profit  P.  =  3  >  P*.   Designate  vertex  4  as  t  . 

Now  let  i  =  2.   Choose  x_  as  indicated.   Assign  t.  to  vertex  8 
which  is  the  endpoint  of  x?.  Using  Step  3  again,  y?  is  selected  to  be 
the  imaginary  branch  ('8,7).   The  gain  g  =  9  and  P  -  P.  +  g  =  12.  Be- 
fore going  further,  check  if  joining  t.  to  t.  would  improve  P*.   With 
(8,1)  as  y  *,  g*  -  9  and  P*  =  P.  +  g  *  -  12.   We  have  to  terminate  the 
,  search  for  x.  's  and  y.  's  since  P*  =  P?  (Step  4).   Had  P*  been  greater 
than  P  we  could  have  continued  by  designating  vertex  7  as  t-  and  look 
for  x^.  and  y„. 

Construct  the  Y's  and  delete  the  X's.   The  network  with  improved 
cost  of  2^9  (not  yet  optimal)  is  shown  in  Fig.  5* 3*  The  cost  is  obtained 
by  subtracting  the  costs  of  x.  and  x?  from  the  initial  value  and  then 
adding  the  costs  of  the  branches  y.  and  y  *. 
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Feasibility  of  the  new  network  is  tested  using  Algorithm  1  since 
we  need  only  to  check  the  branch  connectivity  between  the  vertex-pairs 
1-2  and  4-8.   This  network  is  feasible  so  we  can  iterate  the  procedure 
using  Fig.  5«3  as  the  new  starting  solution. 

Satisfying  the  vertex  requirement  does  not  necessarily  mean  fea- 
sibility is  preserved,   Fig,  5.4  obtained  by  the  optimizing  procedure  is 
not  feasible  as  it  has  a  minimum  connectivity  of  2. 
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V.   CONCLUSION 

Analysis  procedures  provide  convenient  tools  for  evaluating  the 
"strength"  or  "weakness"  of  flow  networks  with  respect  to  some  physically 
meaningful  criteria.   The  computer  analysis  programs  were  prepared  as  im- 
plementation of  the  algorithms  developed  and  are  capable  of  investigating 
large  networks  based  on  either  branch  or  vertex  survivability  measure. 

It  was  shown  that  the  analysis  tools  are  not  only  used  for  measuring 
performance  of  existing  systems  but  are  also  applied  effectively  in  their 
design  particularly  when  either  explicit  solutions  are  impractical  or  so- 
lutions do  not  appear  to  be  obtainable.   A  heuristic  algorithm  for  the 
design  of  minimum-cost  survivable  networks  is  presented.   Analysis  tech- 
niques are  found  indispensable  in  the  procedure. 

Further  study  is  indicated  in  design  perhaps  by  encoding  the  algorithm 
for  computer  running  and  to  investigate  the  following: 

1.  What  level  of  backtracking  would  be  considered  optimal? 

2.  How  often  does  inf easibility  occur  in  the  starting  routine?  The 
probability  of  occurrence  may  be  different  from  that  obtained  by  Steiglitz 
since  the  criteria  are  not  the  same. 

3.  Are  refinements  described  in  [l6l  applicable  to  networks  with 
redundancy  greater  than  2? 
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APPENDIX  A 


FLOWCHARTS 


A-l.   s-t  BRANCH  CONNECTIVITY 


G> 


Read-in  K-matrix, 
source  and  sink 


-yes 


Initialize  flow 
matrix  to  zero 


yes 


Augment  flow  by- 
Augment  ati  on  Rout 


© 


Erase  vertex 
labels 


Initialize  varia- 
bles &  counters 


Search  for  Aug- 
mentation Path  by- 
Labeling  Routine 


no 


no- 


Search  for  br- 
cutset, 
Record  cut  value 


Print  cutset 
and  s-t 

connectivity 
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A-2.      MULTITERMINAL  BRANCH  CONNECTIVITY 
(  Undirected  Network  ) 


C 


& 


XBAR 


nas  single 


List  down  single 
vertex  in  XBAR 
as  processed. 


Read  Connection 
Matrix 


Initialize  vari- 
ables, matrices, 
source,  and  sink 


Find  initial 
s-t  cut. 
Use  Algorithm  1 


List  down  labeled 
vertices  under  X  and 
unlabeled  ones  under 
XBAR. 


Update  Terminal 
Capacity  Matrix 


X  and  XBAR  single 
X  has  single 


X  and  XBAR 
have  >  1 


Condense  XBAR, 
Increment  XBAR- 
counter 


->=. 


List  down  single 
vertex  in  X  as 
processed. 


Choose  source  and 
sink  from  X 


Choose  source  and 
sink  from  XBAR 


«fi£- 


Find  s-t  cut 
Use  Algorithm  1 


0 
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A- 2,  (continued) 


•yes 


Print  Branch- 
connectivity 
matrix 


no- 


Decrement  XBAR 
counter 


Label  XBAR  ver- 
tices, all  other 
vertices  unlabeled 


E> 
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A-3.   MULTITERMINAL  BRANCH  CONNECTIVITY 
(  Directed  Network  ) 


©■ 


© 


-no- 


Choose  a  source 
vertex 


Find  initial 
s-t  cut 


Find  an  Aug- 
mentation Path 
(Routine  A) 


f     Read  N  and 
Connection  Matrix 


Choose  a  sink 
vertex 


—  yes. 


© 


yes- 


Print 
Connectivity 
latrix 


•no' 


Enter  cut  value 
in  Terminal   Cap- 
acity Matrix 


© 
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A-3.   (continued) 


V 


no 


yes 


A  g-t  path  is  found; 
find  a  Homing  Path 
by  Routine  A 


yes 


no 


A  loop  of  flow 
is  detected,  set 
branch  flows  to 
zero. 


Erase  vertex 
labels 


Find  a  Trunca- 
tion Path  by 
Routine  T 


no 


yes 


Trace  and  Erase 


Truncate 


Augment 
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A-4.    s-t  VERTEX  CONNECTIVITY 


©■ 


no 


Initialize  flow 
matrix  to  zero 


no 


Search  for 
vertex  cutset 


Print  s-t 
connectivity 


C 


Read,  Connect- 
ion Matrix 


C 


Read,  Source, 
and  Sink 


Erase  vertex 
labels 


Search  for  Aug- 
mentation Path 
by  Labeling  Rout 


no« 


yes 


yes- 


Augment  flow 
along  this  path 
by  Augmentation 
Routine 


© 
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A-5.      UNIFORM  R-  CONNECTIVITY 


Read,   N,   R, 
connection  matrix 


G> 


Choose  source 
vertex  • 


yes 


Choose  a  sink  vertex 

v..     List  no.    of  ver- 

3 
tices  previously  des- 
ignated sink  that  are 
adjacent  to  v 


Network  is 
R-connected 


no 


Remove  all  the  pre- 
viously designated 
sink  vertices  adja- 
cent to  v .   and  verify 

3 
R-connectivity  be- 
tween v     and  v  .  using 
s  3 

Algorithm  k. 


no 
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APPENDIX   B 


COMPUTER  PROGRAMS 


B-l.   PROGRAM  MAXCON-1 


C    THIS  PROGRAM  CGMPUTES  FOR  THE  MULT  ITERM IN AL   ERANCH- 

C  CONNECTIVITY  OF  A  DIRECTED  OR  MIXED  NETWORK.  THE  NETWORK 

C  IS  REPRESENTED  BY  ITS  CONNECTION  MATRIX  FOR  COMPUTATION. 

C  THE  NETWORK  MUST  HAVE  ONLY  SINGLE  BRANCH  BETWEEN  VERTICES. 

C 

C    THE  PROGRAM  PROCEEDS  BY  READING-IN  THE  INPUT  CATA  CON- 

C  SI  STING  OF  THE  TOTAL  NUMBER  OF  VERTICES  AND  THE  CONNECT- 

C  ION  MATRIX.  IT  THEN  DETERMINES  THE  INITIAL  S'-T  CUT  WITH 

C  THE  FIRST  AND  NTH  VERTICES  AS  SOUkCE  AND  SINK  RESP. 

C  HAVING  FOUND  AN  S'-T  CUT,  IT  THEN  FINDS  OTHER  S-T  CUTS 

C  KEEPING  THE  SINK  VERTEX  CONSTANT  AND  VARYING  ONLY  THE 

C  SOURCE  VERTEX  FROM  2  TO  N.  AFTER  COMPLETING  THIS  LOOP, 

C  THE  SINK  VERTEX  IS  SET  TO  N-l  AND  THEN  FINDS  AN  S'-T  CUT 

C  BETWEEN  VERTICES  1  AND  N-l.  THE  SOURCE  IS  AGAIN  VARIED 

C  FROM  2  TO  N-l  FOR  THE  OTHER  S-T  CUTS, 

C 

C    THE  PROCESS  IS  REPEATED  UNTIL  N(N-l)  CUTS  ARE  DETER- 

C  MINED.  EACH  REPETITION  INVOLVES  SUBROUTINE  FLCVAR. 

C 

C  THE    OUTPUT    IS    A    PRINTOUT    OF    THE    NXN    BRANCH-CONNECTIVITY 

C    MATRIX. 

C 

C IDENTIFIER    NOMENCLATURE    . .     . 

C 

C       ARRAYS 

C  KMAT INPUT    CONNECTION    MATRIX    WITH    ENTRIES    0,     1,    -1 

C  DMAT S'-T    FLOW    MATRIX    REPRESENTING    AN    S'-T    FLOW 

C  PATTERN.    THIS    MATRIX    IS    UPDAThD    BEFORE    A    NEW    SET 

C  OF    S    AND    S'     VERTICES    ARE    DESIGNATED. 

C  FLOMAT-    S-T    FLOW    MATRIX    REPRESENTING    S-T    FLOW    PAT- 

C  TERN    RESULTING    FROM    AUGMENTATION    AND    HOMING    STEPS. 

C  NONZERO    ENTRIES    OF    THIS    MATRIX    UPDATE     CMAT. 

C  TERMAT-    MATRIX    DENOTING    THE    MAXIMUM    BRANCH-CCNNEC T- 

C  IVITY    BETWEEN    ANY    VERTEX-PAIR. 

C 

C       VECTOR 

C  NODE LIST    OF    VERTICES    OF    THE    NETWORK    GRAPH. 

C 

C  ...  ADDITIONAL  INFORMATION  

C    THIS  PROGRAM  IS  CODED  IN  FORTRAN  LANGUAGE.  NETWORK  OF 

C  ANY  SIZE  CAN  BE  HANDLED  BY  THIS  PROGRAM  WITHIN  THE  LIMIT- 

C  ATIONS  OF  THE  COMPUTER  STORAGE.  THE  DIMENSION,  READ,  AND 

C  WRITE  STATEMENTS  MUST  HOWEVER  BE  ADJUSTED  ACCORDINGLY. 

C 

C MAIN  PROGRAM 


C 


C 
C 
C 


INPLICIT  INTEGER(A-Z)  t       r%i  „,oe  or% 

DIMENSION  KMAT(25,25),  FLOMAT( 25 ,25 ) ,  DMAT(25,25) 

DIMENSION  TERMAT ( 25,25) 

DIMENSION  NODE (25) 

COMMON  /ARR1/  KMAT  ,  FLOMAT , DMAT 

COMMON  /ARR2/  TERMAT 

COMMON  /CH/  NODE 

COMMON  /PT/  SOURCE, SPRIME, SINK, AMARK ,BRCUT , N ,G 
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READ  (5,9)  N 

READ  (5,10)  ( (KM4T(I,J)  ,J=1  ,25)  ,I=1,N) 
C  INITIALIZE  TERMINAL  MATRIX 

C 

DO  2  X=1,N 
C 

DC  1  Y=l,25 

TERMAT( X,Y)=0 

1  CONTINUE 
C 

2  CONTINUE 
C 

C 

C  COMPUTE  FOR  THE  N(N-l)  BRANCH  CUTS 

C 

DC  8  K=1,N 

J=N+1-K 
C 

C  SET    TO    ZERO    THE    FLOW    MATRICES 

C 

DO   4    L=1,N 
C 

DO    3   M=1.N 

DMAT(L,M)-0 

FLCMAT(L,M)=0 

3  CONTINUE 


C 

c 
c 


4    CONTINUE 
SPRIME=0 

DO    7    1=1, N 

IF     (I  .EQ.J)    GC   TO    7 

SOURCE = I 

SINK=J 

BRCUT=0 

CALL    FLCVAR 
C 
C  ENTER    CUT    VALUE    INTO    THE    TERMINAL    MATRIX 

TERMATU  ,J)=BRCUT 

SPRIME=1 
C 
C  TRANSFER    FLOMAT    ENTRIES     INTO    DMAT 


C 


C 
C 
C 


DO    6    L=1,N 


DO    5    M=1,N 

IF    (FLCMAT(L,M) .EQ.O)     GO    TO    5 

DMAT(L,M) =FLOMAT(L,M) 
C 
C  RESET    FLOMAT    TO    ZERO 

FLOMAK  L,M)=0 
5    CONTINUE 


6  CONTINUE 

7  CONTINUE 

8  CONTINUE 


C 
C 

C  COMPUTATION  COMPLETED 

WRITE  (6,11) 

WRITE  (6,12) 

CALL  PRTMAT 

STCP 
C 

9  FORMAT  (16) 

10  FORMAT  (2513) 

11  FCRMAT  (•Ot ,T14,» BRANCH-CONNECTIVITY  MATRIX*) 

12  FORMAT  (T14,' ') 

END 
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C  SUBROUTINE  FLOVAR   

C 

C    THIS  SUBROUTINE  DETERMINES  AN  S-T  CUT  ONCE  AN  INITIAL 

C  S«-T  CUT  IS  FOUND.  IT  MAXIMIZES  AN  S-T  FLOW  USING  THE 

C  EXISTING  S«-T  FLOW  PATTERN. 

C 

SUBROUTINE  FLOVAR 

IMPLICIT  INTEGER(A-Z) 

DIMENSION  KMAT(25,25),  FLOMAT ( 25  ,25) ,  DMAT(25,25) 

DIMENSION  NJDE(25) 

CCMMQN  /ARR1/  KM AT , FLOM AT t DMAT 

COMMON  /CH/  NODE 

COMMON  /PT/  SOURCE, SPRIME, SINK, AMARK, BRCUT , N,G 

c 

C  ERASE  ALL  VERTEX  LABELS 

C 

1  DC  2  J  =  1,N 
NQDEt JJ=0 

2  CONTINUE 
C 

C 

C  SEARCH  FOR  AUGMENTATION  PATH 

CALL  AUGMNT 

IF  (AMARK)  3,24,15 

C 

C  AN  S-G  AUGMENTATION  PATH  IS  FCUND 

C  SEARCH  FOR  HOMING  PATH 

3  Q=G 
C 

A  DO  5  J=1,N 

IF  (DMAT(Q,J).EQ.l)  GO  TO  6 

5  CONTINUE 
C 

GO  TO  1 

6  NODE(J)=Q 
Q  =  J 

IF  (Q.EG.G)  GO  TO  19 

IF  (J. EU. SINK)  GO  TO  7 

GO  TO  4 
C 

C  A  HCMING  PATH  IS  FOUND 

C  SEARCH  FOR  TRUNCATION  PATH 

7  Z  =  G 

IF    (G.EQ. SPRIME)     GO    TO    12 

C 

8  DO    9    1=1, N 

IF    (DMATd  ,Z).EQ.l)    GO    TO    10 

9  CONTINUE 
C 

IF    (Z.EG.G)    GO   TO    12 
GO    TO    11 

10  NCDE(I)=Z 

IF    (Z.EQ.G)    GO    TO    21 

IF    (Z.EQ.SPRIME)     GO    TO    11 

GO    TO    8 

C  A    TRUNCATION    PATH    IS    FOUND 

C  TRUNCATE    SPRIME-G    PATH 

11  I=NODE(Z) 
DMATIZ, I)=0 
Z=I 

IF    (Z.EQ.G)    GO    TO    12 
GC    TO    11 

C  TRACE    AND    ERASE    G-T    PATH 

12  Q=SINK 

13  J=NODE( G) 
FLOMAT(JiQ)=l 
OMAT( J,Q)=0 

Q  =  J 

IF    (J.EQ.G)    GO    TO    14 
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GO  TO  13 
C 
C  AUGMENT  S-G  PATH 

14  IF  (G.EQ. SOURCE)  GC  TO  23 
J  =  G 

GO  TO  16 
C 
C  AUGMENT  S-T  PATH 

15  J=SINK 

16  I=NODE(J) 

IF  (I .LT.O)  GO  TO  17 
C 
C  INCREASE  FORWARD  FLOW  BY  UNITY 

FLOMAT<  It  J J=l 

GC  TO  18 
C 
C  DECREASE  BACKFLOW  BY  UNITY 

17  I=~I 
FLCMAT(  J, I  )=0 

18  J=I 

IF  (J.EQ.SCURCE)  GO  TO  23 

GO  TO  16 
C 
C  A  LOOP  OF  FLOW  IS  FOUND  BY  HOMING  ROUTINE 

19  Q  =  G 

20  J=NODE(G) 
DMAT( J,Q)=0 
Q  =  J 

IF  (Q.EQ.G)  GO  TO  1 

GO  TO  2  0 
C 
C  A  LOOP  OF  FLOW  IS  FOUND  BY  TRUNCATION  ROUTINE 

21  Z=G 

22  I=NODE(Z) 
DMATlZt  I)=0 
Z=I 

IF  (Z.EQ.G)  GO  TO  1 
GC  TO  22 

23  BRCUT=BRCUT+1 
GO  TO  1 

C 

C  NO  MORE  AUGMENTATION  PATH  IS  FOUND.  FLCW  PATTERN 

C  IS  MAX 

24  RETURN 
END 

C 
C 

c 

C  SUBROUTINE  AUGMNT   . 

C 

C  THIS    SUBROUTINE    IS    USED    TO    DETECT    AN    S-G    OR    S-T    PATH 

C    USING    THE    MODIFIED    FORD    AND    FULKERSON    LABELING    ALGORITHM 

C 

SUBROUTINE  AUGMNT 

IMPLICIT  INTEGER(A-Z) 

DIMENSION    KMAT(25,25i,    FLOMAT ( 25 , 25 ) ,    0MAT(25,25) 

DIMENSION    N0DE(25i,    JST0R5U0),     ISTORE(IO) 

CCMMON    /ARR1/    KMAT , FLOMAT , DMAT 

CCMMON    /CH/    NODE 

CCMMON    /PT/    SOURCE, SPRI ME t SINK t AMARK , BRCUT  ,  N, G 
C 

I =SOURCE 

L=0 

M=0 

Z=0 

AMARK=0 

NCDE( I)=SOURCE      . 
C 
C  SEARCH    FOR    AUGMENTING    PATH 


C 


DO    4    J=1,N 
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C      tc  / JS?I, FCR  CONNECTING  BRANCH  BETWEEN  VERTEX  PAIRS 

IF  (DMAT( J,I J.EQ.OJ  GO  TO  2 

IF  (J .EQ. SOURCE)  GO  TO  4 
C 

C      TC  /^f^1!1^  BACKFLOW  IS  ENCOUNTERED 

IF  (NODt(JJ.NE.O)  GO  TO  4 
C 
C  LABEL  VERTEX  J 

NCDE( J)=-I 

GO  TO  3 

2  IF  (KMAT(-ItJ).NE.l)  GO  TO  4 
C 

C  A  FORWARD  BRANCH  IS  ENCOUNTERED.  TEST  FCR  FORWARD 

C  FLOW 

IF  <DMATU,J)  .NE.O)  GO  TO  7 

IF  (FLOMAT(l,JJ.NE.O)  GO  TO  4 

IF  (NODE(J).NE.O)  GO  TO  4 
C 
C  LABEL  VERTEX  J 

NODE(J)=I 
C 
C  IS  THE  SINK  VERTEX  REACHED? 

3  IF  (J.EQ.SINK)  GO  TO  8 
M  =  M+1 

JSTORE(M)=J 

4  CONTINUE 
C 

IF  (Z.NE.O)  GO  TO  6 

IF  (M.EQ.O)  GO  TO  9 
C 

C  TRANSFER  CONTENTS  OF  JSTORE  TO  ISTORE 


C 


DC  5  S  =  1,M 

ISTORE( S)=JSTGRE( S) 

CONTINUE 


Z  =  M 

M  =  0 

6  I=ISTORE(Z) 
Z=Z-1 

GO  TO  1 
C 
C  AN  S-G  AUGMENTATION  PATH  IS  FCUND 

7  G=l 
AKARKs-1 
GO  TO  10 

C 

C  AN    S-T    AUGMENTATION    PATH    IS     FOUND 

8  AMARK=1 
GO    TO    LO 

C 

C  NO  AUGMENTATION  PATH  FOUND 

9  AMARK=0 
10  RETURN 

END 
C 
C 

c 

C  SUBROUTINE  PRTMAT   

C 

C      THIS  SUBROUTINE  PRINTS  TERMAT  IN  PROPER  FCRM  AND  IN 

C  FCRMAT  3513. 

C 

C    IF  ARRAY  LENGTH  IS  .LE.  35  A  SINGLE  LOOP  SECTION  PRINTS 

C  THE  MATRIX,  OTHERWISE  IT  PRINTS  THE  FIRST  35  CCLUMNS  FGL- 

C  LCWED  BY  THE  NEXT  35  COLUMNS  AND  SO  ON  IN  MULTIPLES  OF  35 

C  USING  MULTIPLE  LOOPS.  THE  REMAINDER  COLUMNS  ARE  PRINTED 

C  IMMEDIATELY  BELOW  THE  LAST  35-COLUMN  GROUP. 

C 

c 
c 
c 
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SUBROUTINE    PRTMAT 

IMPLICIT     INTEGER(A-Z) 

DIMENSION    TERMAT(25,25) 

CCMMON    /ARR2/    TERMAT 

CCMMCN    /PT/    SOURCE, SPRIME,SINK, AMARK, BRCUT , N,G 

C  HOW    MANY    MULTIPLES    OF    3  5    IS    N? 

MULPLE=N/35 

IF    (MULPLE.GE.l)     GO    TO    2 
C 
C  THIS    SECTION    PRINTS    ARRAY    WITH    LENGTH    .LT.    35 


C 


C 
C 


DO    1    1=1, N 

WRITE     (6,8)     I, (TERMAT(I,J), J=1,N) 
1    CONTINUE 


GO  TO  7 
C 

C  THIS  SECTION  PRINTS  ARRAY  WITH  LENGTH  • GE .  35 

C  IN  GROUPS  OF  35 

2  M  =  l 
K=l 

3  MPLS35=M+35 
C 

00    4    I=1,N 

WRITE     (6,8)     I, (TERMAT(I,J J, J=M,MPLS35) 

4  CONTINUE 
C 

M=M+35 

IF     (K.EQ.MULPLE)     GO    TO    5 

WRITE     (6,9)     K 

WRITE     (6,10) 

K  =  K+ 1 

GO    TO    3 
C 
C  THIS    SECTION    PRINTS    REMAINDER    COLUMNS 

5  REMDR=N-MULPLE*3  5 

IF  (REMDR.LE.O)  GO  TO  7 
WRITE  (6,11) 
WRITE  (6,12) 
C 

DO  6  L=1,N 

WRITE  (6,8)  L, (TERMAT(LtJ) , J=M,N) 

6  CONTINUE 


7  RETURN 

8  FORMAT  CO',  T8,  12,  •/',  IX,  3513) 

9  FCRMAT  (//T14,« CONTINUATION ',  13,  2X,  *0F  TERMAT*) 
1C  FORMAT  (714, « •) 

11  FORMAT  (//T14,'LAST  SECTION  UF  TERMAT1) 

12  FCRMAT  (T14,» * «) 

END 
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B-2.   PROGRAM  MAXCON-1A 


C    THIS  PROGRAM  COMPUTES  FCR  THE  BRANCH- CONNECT  I VITY  BET- 

C  WEEN  ANY  TWO  GIVEN  VERTICES  OF  A  DIRECTED  OR  MIXED  NET- 

C  WCRK.  THE  NETWORK  GRAPH  IS  REPRESENTED  BY  ITS  CONNECTION 

C  MATRIX  FOR  COMPUTATION.  NETWORK  MUST  HAVE  ONLY  SINGLE 

C  BRANCH  BETWEEN  VERTICES. 

C 

C    THE  NUMBER  CF  VERTICES,  THE  SELECTEC  SOURCE  ANC  SINK 

C  VERTICES,  AND  THE  CONNECTION  MATRIX  ARE  THE  INPUTS.  THE 

C  PKCGKAM  DOES  A  MAXFLOW  CALCULATION  USING  THE  MOCIFIED 

C  FORD  AND  FULKERSON  LABELING  ALGORITHM.  COMPUTED  CUT  VALUE 

C  IS  THE  NETWORK  CONNECTIVITY  WITH  RESPECT  TO  THE  SPECIFIED 

C  VERTEX-PAIR. 

C 

C         SEVERAL    VERTEX-PAIRS    MAYBE    LISTED    WITH    THE     INPUT    AND 

C    THE     PROGRAM    DOES    SIMILAR    MAXFLOW    CALCULATION     FOR    EACH 

C    PAIR. 

C 

C    THE  OUTPUT  IS  A  PRINTOUT  CONSISTING  OF  THE  GIVEN  SOURCE 

C  AND  SINK,  THE  CONNECTIVITY  VALUE,  AND  THE  BRANCH-CUTSET 

C  AS  DENOTED  BY  THE  VERTICES  TO  WHICH  EACH  BRANCH  CUT  IS 

C  CONNECTED. 

C 

C    THE  PRESENT  DIMENSION  AND  READ  STATEMENTS  MAYBE  ADJUST- 

C  ED  ACCORDING  TO  THE  SIZE  OF  THE  NETWORK  TO  BE  ANALYZED. 

C 

C  THIS    PROGRAM    MAYBE    PREFERRED    OVER    MAXCON-1    WHEN    ONLY    A 

C    FEW    SOURCE-SINK    PAIRS    ARE    TO    BE    ANALYZED. 

C 

c 

IMPLICIT     INTEGER(A-Z) 

DIMENSION    KMAT(50,50)t    FLOMAT (50  ,50)  t    N0DE(50),     NX(5G) 
DIMENSION    N0TNX(20),     JSTORE(203,     IST0RE(2C) 
C 

READ    (5,22)    N 

READ    (5,23)     (  (KM  AT ( I , J ) , J= 1 , 50 ) , I  =  1 , N ) 

WRITE     (6,25) 

1  READ    (5,24)     SOURCE,  SINK 
IF     (SOURCE. EQ.O)     GO    TO    21 

C 

C  INITIALIZE    FLOW    MATRIX    TO    ZERO 

C 

DO    3    I=1,N 
C 

DC    2    J=1,N    ' 

FLOMAT( I , J)=0 

2  CONTINUE 
C 

3  CONTINUE 
C 

c 

C  ERASE    VERTEX    LABELS 


C 


C 
C 


c 


4  DC    5    J=1,N 
NGDE( J) =0 

5  CONTINUE 


I=SOURCE 
L=0 

M  =  0 

z  =  c 

BRCUT=0 

NCDE( I) =SOURCE      - 

C  SEARCH    FOR    AUGMENTING    PATH 


DO    10    J=1,N 
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C      r,     ,J£SJ.frCR    CONNECTING  BRANCH  BETWEEN  VEPTEX-PAIR 

IF     (KMAM  I,  J  )  )     7,  10,8 
C 

C  A    BACKWARD    BRANCH    IS    ENCOUNTERED.    TEST     FCR    PCSI- 

C  TIVE    6ACKFL0W. 

7  IF     (FL0MAT1 J, I) .ECO)     GO    TO    LO 
IF    (NODEUJ.NE.O)     GO    TO    10 

C 

C  LABEL    NODE(J) 

NODE( J)  =-I 

GO    TO    9 
C 

C  A    FGRWARD    BRANCH    IS    ENCOUNTERED.    TEST    FCR    ZERO 

C  FORWARD    FLOW. 

8  IF     (FLOMAT (I, J) .NE.O)     GO    TO     10 
C 

C  IS    NODE(J)    LABELED? 

IF     (NQDE(J) .NE.O)     GO    TO    10 

NCDE{ Ji  =1 
C 
C  IS    THIS    THE    SINK    VFRTEX? 

9  IF    (J.EQ.SINK)    GO    TO    13 
M=M+1 

JSTORE(M)=J 
10    CONTINUE 
C 

IF     (Z.NF.O)     GO    TO    12 

IF     (M.EQ.O)     GO    TO    16 
C 
C  TRANSFER    CONTENTS    OF    JSTCRE    TO    ISTORE 


C 


CO  11  S  =  1,M 
ISTOP.E(S)=JSTCRt(S) 
11  CONTINUE 


Z=M 
M  =  0 

12  I=ISTOPf-£Z) 
Z=Z-1 

GC  TO  6 
C 

C  WE  REACH  HERE  IF  AN  AUGMENTING  PATH  IS  FOUND 

C  ADJUST  FLOW  ON  ALL  BRANCHES  ALONG  TFIS  PATH, 

13  I=NODE(J) 

IF    (I .LT.O)     GO    TO    14 
C 
C  INCREASE    FORWARD    FLOW    BY    UNITY 

FLCMAT(  I ? Ji  =1 

GO  TO  15 
C 
C  DECREASE  BACKFLCW  BY  UNITY 

14  I=-I 

FLOMAT (J,  I)  =  0 

15  J  =  I 

IF  (J. EO. SOURCE)  GO  TO  4 

GC  TO  13 
C 

C  THIS  STEP  IS  REACHED  IF  LABELING  AND  FLCW  AUGMENT- 

C  ATION  IS  COMPLETED.  SEARCH  FOR  BRANCH  CUTS. 

C 

16  DC  19  I=1,N 

IF  (NODEU  J.EQ.O)  GO  TO  19 
C 

DO  18  J=1,N 

ACD  =  1 

IF  (KMAT(  I  tJKEQ.O)  GC  TO  18 

IF  (NODE( J) .NE.Oi  GO  TO  18 
C 
C  IF  FLOW  CCCURS  IN  BOTH  DIRECTIONS,  NET  FLCW  IS  ZERO 

IF  (FLOMATd, J)  .EQ.FLOMAT(  J, I)  )  GO  TO  18 

IF  (FLOMATd, J)  .EQ.li  GO  TO  17 

ADO=FLOMAT( I , J) 
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17  BRCUT=BRCUT+ADD 
L=L+1 

NX(L)=I 
NCTNX (L)=J 

18  CONTINUE 
C 

19  CONTINUE 
C 

c 

C  COMPUTATION  TERMINATES  SINCE  NO  MORE  AUGMENTING 

C  PATH  IS  FOUND. 

IF  (BRCUT.NE.O)  GO  TO  20 

WRITE  (6,26)  SOURCE, SINK, BRCUT 

GO  TO  1 

20  WRITE  (6,27)  SOURCE , S INK,  BRCUT ,  < NX( J  J ,NOTNX ( J ) , J  =  l , L 
GO  TO  1 

21  STOP 
C 

22  FORMAT  (16) 

23  FORMAT  (2513, /25I3) 

24  FCRMAT  (2  16) 

25  FORMAT     { « C ,     10X,     ' SOURCE* ,     10X,     'SINK',     10X,     '  BR- 
1  CONNECTIVITY',     10XT     '8R-CUTSETM 

26  FORMAT     ({0«,     12X,     12,     13X  *     12,     17X,     12) 

27  FORMAT    CD1,     12X,    12,     13X,    12,    17X,    12, 
1  (T67,     13,     IX,      »~ ',     13/)) 

END 
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B-3.   PROGRAM  MAXCON-2 


C     THIS  PROGRAM  COMPUTES  FOR  THE  MU LT  ITER  MI NAL  BRANCH- 

C  CONNECTIVITY  OF  ANY  GIVEN  UNDIRECTED  NETWORK.  TFE  NETWORK 

C  GRAPH  IS  PRESENTED  FOR  COMPUTATION  BY  ITS  CONNECTION  MAT- 

C  RIX.  NETWORK  CAN  HAVE  PARALLEL  BRANCHES  BETWEEN  VERTICES. 

C 

C     THE  PROGRAM  PROCEEDS  BY  FIRST  PEADING-IN  THE  NUMBER  OF 

C  VERTICES,  THE  INITIAL  VERTEX-PAIR  CHOSEN  ARBITRARILY, 

C  ANC  THE  CONNECTION  MATRIX.  THE  CUT  BETWEEN  THE  TWO 

C  SELECTED  VERTICES  IS  OBTAINED  USING  SUBROUTINE  MINCUT. 

C  THE  CONNECTIVITY  FOR  THE  CUT(X,XBAR)  IS  THE  COMPUTED  CUT 

C  VALUE.  THIS  IS  ENTERED  INTO  THE  TERMINAL  CAPACITY  MATRIX 

C  BY  SUBROUTINE  MATRXT.  THE  SUCCEEDING  VERTEX-PAIRS  ARE 

C  SELECTED  BY  THE  PROGRAM  FROM  THE  GENERALIZED  X-VERTEX- 

C  GROUP.  SUBROUTINE  MINCUT  IS  EMPLOYED  EACH  TIME  FOR 

C  OBTAINING  THE  CUT  VALUE.  WHEN  ALL  VERTICES  IN  X  ARE  EX- 

C  HAUSTED,  THE  X-VE RTE X-GROUP  IS  CONDENSED  INTC  A  SINGLE 

C  VERTEX  BY  SUBROUTINE  MATRXC  AND  THE  X6AF- VER T E X-GROUP  IS 

C  IN  TURN  PROCESSED  SIMILARLY.  FOR  EACH.  COMPUTATION,  SUB- 

C  ROUTINE  MATRXT  UPDATES  THE  TERMINAL  CAPACITY  MATRIX. 

C 

C    V.ITH  N  AS  THE  TOTAL  NUMBER  OF  VERTICES,  THE  PROGRAM 

C  WILL  CO  N-l  SUCH  COMPUTATIONS  AND  FINALLY  PRESENT  A  MULTI- 

C  TERMINAL  BR ANCH-CCNN ECT IV  IT Y  MATRIX  AS  OUTPUT. 

C 

C  IDENTIFIER  NOMENCLATURE  

C 

C       ARRAYS 

C  KMAT INPUT    CONNECTION    MATRIX    WHICH    HAS    ENTRIES 

C  GREATER    THAN    UNITY     INDICATING    PARALLEL    BRAN- 

C  CHES    BETWEEN    VERTICES 

C  CMAT-^ — MATRIX    REPRESENTING    THE    NEW    GRAPH    FORMED    RES- 

C  ULT1NG    FROM    VERTEX    CONDENSATION.     IT    HAS    ALSO 

C  ENTRIES    GREATER    THAN    UNITY. 

C  FLOMAT--MATRIX    INDICATING    MAXIMUM    FLOW    THAT    CAN    EXIST 

C  BETWEEN    ANY    VERTEX-PAIR    ALLOWING    UNIT    CAPA- 

C  CITY    FOR    EACH    BRANCH. 

C  TERMAT — MATRIX    DENOTING    THE    MAXIMUM    CONNECTIVITY    BET- 

C  WEEN  ANY  VERTEX-PAIR. 

C      NCTVTXJ -ARRAY  STORAGE  FOR  SEQUENCE  OF  CUTS (X, XBAR ) * 
C  HAVING  MORE  THAN  ONE  VERTEX  ON  EACH  SIDE. 

C 

c  NX^ --VECTOR  STORAGE  FOR  X-VERTICES  IN  CUT(X,XBAR) 

C  CURRENTLY  BEING  PROCESSED. 

C  NGTNX VECTOR  STORAGE  FOR  XBA R-VERT ICES  IN  CUT 

C  (X,XBAR)  CURRENTLY  BEING  PROCESSED. 

C  NLIST— -STORAGE  FOR  VERTICES  WHICH  HAVE  OCCURED  SING- 

C  ULAR  IN  A  CUT{X,XBAR). 

C 

C   CONSTANTS  AND  COUNTERS 

c  N TOTAL  NUMBER  OF  VERTICES  IN  THE  GRAPH 

C  SOURCE — CHOSEN  SOURCE  VERTEX 

C  SINK CHCSEN  SINK  VERTEX 

£  RCW DESIGNATES  ROW  NUMBER  IN  ARRAYS  VRTX  AND 

C  NOTVTX 

C  BRCUT DENOTES  CUT  VALUE 

C  CCUNT1— NUMBER  OF  VERTICES  IN  NX 

C  CCUNT2 — NUMBER  OF  VERTICES  IN  NCTNX 

C  NCOUNT — RUNNING  TOTAL  OF  VERTICES  IN  NLIST 

C  ADDITIONAL  INFORMATION  

C     THIS  PROGRAM  IS  CODED  IN  FORTRAN  LANGUAGE.  THE  PRESENT 

C  DIMENSION  STATEMENTS  REQUIRE  I50K  OF  COMPUTER  STORAGE. 

C  LARGER  NETWORKS  CAN  BE  HANDLED  BY  ADJUSTING  ThE  DIMENSION, 

C  REAC,  AND  WRITE  STATEMENTS. 

C 
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C .  .  .  .   MAIN  PROGRAM 

C 


c 
c 


CCMMON  /ARRl./  KMAT,f 
CCCMON  /ARR2/  TERMAI 
CCMMON    /Chi/    NODE,NLIST 


CCMMON    /CH2/    NX,NOTNX 

CCMMON    /PT1/     N, SOURCE, SINK, ROW, MARK, FLAG, BRCUT 

CCMMON     /PI 2/    COUNT1  ,COUNT2 , NCCUNT 

PEAD    (5,29)     N, SOURCE, SINK 

READ    (5,30)     (  (KMAT(  I  ,J)  ,J  =  1  ,35.)  ,I=1,N) 

C 

C  INITIALIZE  VARIABLES 


DO  3  I =1,N 

DC  2  K=l,10 

DO    1    J=1,N 
CMAT( I , J)=KMAT( I, J  J 
VRTXd  ,  J)  =0 
NOTVTX( I,J)=0 
NLIST ( J)=0 

1  CONTINUE 

2  CCKTINUE 

3  CONTINUE 

DO    5    1=1,35 

DO    4    J  =  l  ,35 
TERMAT(  I  ,J)=0 

4  CONTINUE 

5  CONTINUE 


BRCUT=0 

NCCUNT=0 

FLAG=0 

MAPK=0 

ROW  =  0 
C 
C  FIND    INITIAL     S-T    CUT 

6  CALL  MINCUT 

C  LIST  DOWN  LABELED  VERTICES  UNDER  NX  AND  UNLABELED 

C  ONES  UNDER  NCTNX 

CCUNT1=0 

CCUNT2=0 
C 

DO  7  J=1,N 

NX(J)=0 

NCTNXt  J)=0 

7  CCNTINUE 
C 

c 

DO  14  I  =  1,N 

C  TEST  IF  THIS  VERTEX  IS  ALREADY  IN  NLIST 

IF  (NCOUNT.EQ.O)  GO  TO  9 
C 

DO  8  JJ=1,NC0UNT 

IF  (NLISTt JJJ .EQ.I)  GO  TO  14 

8  CONTINUE 
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c 
c 

C  TEST    IF    IT    BELONGS    TO    SOME    CONDENSED    NCTVTX 

9     ISTCRE=Rr 
1C    IF    (ISTORE.EQ.O)     GC    TO    12 


C 

c 


c 


c 

c 


c 


DC    11    KK=ltN 

IF    (NOTVTX(ISTORE,KK) . EQ. I)     GO    TC    14 
11    CONTINUE 


I STORE=ISTOkE-1 
GC    TO    10 

12  IF     <NODE(  D.EQ.O)     GO    TO    13 
CGUNT1=C0UNT1+1 
NX(C0UNT1  )=I 

GC    TO    14 

13  CCUNT2=C0UNT2+1 
NC7NX(C0UNT2)=I 

14  CONTINUE 
C 

c 

C  UPDATE  TERMINAL  CAPACITY  MATRIX  ENTRIES 

CALL  MATRXT 
C 

IF  (C0UNT1.FQ.1)  GC  TO  18 

IF  (C0UNT2.EQ.1)  GO  TO  17 
C 
C  BOTH  NX  AND  NOTNX  HAVE  MORE  THAN  ONE  VERTEX 

MARK---0 

RCW  =  RC1W  +  1 
C 
C  STORE    NX    AND    NOTNX    IN    VRTX    AND    NOTVTX    RESP 


DC    15    LL=1,C0UNT1 
VRTX(ROW,LL)=NX(LL  ) 
15    CCNTINUE 


DC    16    LL=1,C0UNT2 

NGTVTX(  ROW  ,LL)  =NCTi\X(LL) 

16  CONTINUE 
C 

c 

C  CONDENSE    NOTVTX     INTO    SINGLE    NODE 

CALL    MATRXC 
C 
C  SELECT    TWO   VERTICES     IN    NX    AS    SOURCE    AND    SINK 

SOURCE=NX( 1) 

SINK=NX(C0UNT1) 

GC    TO    6 
C 
C  NX    HAS    MORE    THAN    ONE    VERTEX    AND    NOTNX    HAS    ONE 

17  NCCUNT=NCCUNT+1 
NLIST(  NCCUNT)=NOTNXm 
SCUPCE=NX(1) 
SINK=NX(CGUNT1) 

GC    TO    6 
r 

C  NX  HAS  SINGLE  VERTEX  AND  NOTNX  HAS  MCPE  THAN  ONE 

18  NCCUNT=NC0UNT+1 
NLIST (NCCUNT)=NX(1 ) 
IF  (C0UNT2.EQ.1J  GO  TO  19 
SCURCE=NGTNX( 1) 
SINK=N0TNX(C0LNT2) 
GO  TO  6 

C 

C  BOTH    NX    AND    NOTNX    HAVE    SINGLE    VERTEX 

19  DC    20    LL=1»NC0UNT  ^^    ^n 

IF     (NLIST(LL ).EQ.NOTNX( 1) )     GO    TO    ^1 
2C    CCNTINUE 
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NCCUNT=NC0UNT+1 

NLI ST( NCCUNT) =NOTNX (1 ) 

C  TEST  IF  ALL  VERTICES  ARE  IN  NLIST 

21  IF  (NCOUNT.EQ.N)  GO  TO  28 

r  TSJ?w?£CII0N  PROCESSES  THOSE  NOTVTX  GROUPS 

C  PREVIOUSLY  CONDENSED 

CCUNT1=0 

CCUNT2=0 
C 

CO  23  MM=1,N 

IF  (NOTVTX(ROW,MM) .ECUO)  GC  TO  22 

COUNT 2= C CUNT  2+1 

22  IF  (VRTX(RCWtMM) .EQ.O)  GO  TO  23 
CCUNT1=CCUNT1+1 

23  CONTINUE 
C 

C 

C  CONDENSE  VRTX  INTO  SINGLE  VERTEX 

C  RELABEL    VERiICES    ACCORDING    TO    PRESENT    CUT 


C 


c 
c 


c 

c 


DO  24  J  =  1,N 
NCDE( J)=0 
24  CONTINUE 


DC  26  I=lfCOUNT2 

DO  25  J=ltN 

IF  (J.NE.NOTVTX(ROW,I J  J  GO  TO  2  5 

NODc( JJ=1 

GC  TO  2o 

25  CONTINUE 

26  CONTINUE 


SOURCE=NOTVTX(ROW,  1) 

S  I  NK  =  NOT VTX  {  RU  ,  COUNT  2  ) 
C 

C  ERASE  PRESENT  VRTX  AND  NOTVTX  ROWS 

C 

DO  27  KK=1  ,N 

VRTXCROWtKK)=0 

NCTVTX(ROW,KK)=0 

27  CONTINUE 
C 

MARK=1 

FLAG=1 

RCW=ROW-l 

CALL  MATRXC 

GO  TO  6 
C 
C  COMPUTATION  TERMINATES 

28  WRITE  { 6,31) 
WRITE  (6,32) 
CALL  PRTMAT 
STOP 

C 

29  FCRMAT  (316) 
3C  FORMAT  (3512) 

31  FORMAT  ( «0' ,T14, 'BRANCH-CONNECTIVITY  MATRIX') 

32  FCRMAT  (T14,» «) 

END 

C 
C 

C  SUBROUTINE  MINCUT  

C 

C     THIS  SUBROUTINE  COMPUTES  FOR  THE  MAXIMUM  CONNECTIVITY 

C  BETWEEN  TWO  SELECTED  VERTEX-PAIR  USING  THE  MCCIFIED  FCRC 

C  AND  FULKERSON  LABELING  ALGORITHM. 

C 
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c 

c. 
c 


c 
c 
c 

c 

c 


c 

c 


c 
c 


c 
c 
c 


c 
c 
c 


FL0MAT(35,35) t 

N3TVTX(35,25) , 

IST0RE(2Q) 


SUBROUTINE    MINCUT 
IMPLICIT     INTEGER(A-Z) 
DIMENSION    KMAT(35 ,35) , 
CIMENSION    VRTX(35t  35) , 
D1MENSI CN    J  ST  ORE (20)  , 
DIMENSION    NLIST( 35) 

/ARR1/    KM AT, FLOMAT,CMAT, VRTX,NOTVTX 

/CHI/    NODE.NLIST 

/PT1/ 

/PT2/ 

INITIALIZE  FLOW  MATRIX 


CMAT(35,35) 

NODE (35) 


CCKMGN 
COMMON 
CCMMON 
CCMMON 


'  KMAT, FLOMAT,CMAT, VRTX,NOTVTX 
NODE,NLI ST 

N, SOURCE , SINK, ROW, MARK, FLAG, BRCUT 
COUNT  1, COUNT 2, NCOU NT 


DO  2  1=1, N 

DC  1  J=1,N 
FLOMAT( I, J)=0 

1  CONTINUE 

2  CONTINUE 


c 

c 

c. 

ERASE    VERT! 

c 

3 

DC    4    J=1,N 
NODEl J) =0 

A 

CONT  INUE 

c 

c 

I=SOURCE 
M=0 

z=o 

BRCUT=0 

NCCE(  I)=SOURCE 

EX  LABELS 


SEARCH  FOR  AUGMENTING  PATH 

5  DO  6  J=1,N 

TEST  FOR  CONNECTING  BRANCH  BETWEEN  VERTEX-PAIR 
IF  (CMATU  ,J)-EO.O)  GO  TO  6 
IF  (FLOMAT(I,J) .EQ.CMATt I, J) )  GO  TO  6 

IS  NODE(J)  LABELED? 
IF  (NODE(J).NE.O)  GO  TO  6 
NCDE( J)=I 

IS  THIS  THE  SINK  VERTEX? 
IF  ( J.EQ.SINK)  GO  TO  9 
M=M+1 
JSTCRE(M)=J 

6  CONTINUE 

IF  (Z.NE.O)  GO  TO  8 
IF  (M.EQ.O)  GO  TO  LO 

TRANSFER  CONTENTS  OF  JSTCRE  TO  ISTORE 

DO  7  S=1,M 
ISTORE(S)=JSTCRE(S) 

7  CONTINUE 

Z  =  M 
M  =  0 

8  I=ISTORE(Z) 
Z=Z-1 

GO  TO  5 

WE  RFACH  HERE  IF  AN  AUGMENTING  PATH  IS  FOUND 
ADJUST  FLOW  ON  ALL  BRANCHES  ALONG  THIS  PATH. 

9  I=NODE(J)  T   ,.  , 
FLOMAT( I,J)=FL0MAT(I,J)+1 
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c 
c. 
c 
c 


J  =  I 

IF     (J 

GO    TO 


EQ.SOURCF) 
9 


GO    TO    3 


THIS  STFP  IS  REACHED  TF  LABELING 
ATION  IS  CnMPLETED.  SEARCH  FOR 


AND  FLOW  AUGMENT- 
BRANCH  CUTS . 


10 


DO 
I  F 


13  I  =1  ,  N 
(MODE(  IJ.EQ.O) 


J  =1 ,  N 


GO  TO  13 


C 

c 


c 
c 


c 
c 
c 

c 
c 
c 
c 
c 
c 


DO  12 

ADD  =  1 

IF  (CMAT(I,J ).E9.0)  GO  T3  12 

IF  (NODE( J) .NF.O)  GT  t^  12 


IF  FLOW  OCCURS  IN  BOTH 


IF  (FLOMATd  »  J) 
IF  (FLOMATd,  J) 

ADD=FLOMAT (T,J| 
BRCIJT=pjrcut+ADD 

CONTINUE 


EQ.FLOMAT( J,I )  ) 
EQ.l)     GO    TO    11* 


DIRECTIONS,     MET 
GO    TO    12 


FLOW    IS    ZERO 


13    CONTINUE 

RETURN 
FND 


SUBROUTINE    MATRXT 


THIS    SUBROUTINE    ENTERS    COMMUTED    CUT    VALUE     INTO    &     TFRMi- 
N6L    CAPACITY    MATRIX     tcr^M     £N0     FURTHER    UPD*TCS     OTHER     EN- 
TRIES   R^SFD    nH    MAXIMUM    Fl  nw    ALONG    A     GIVEN!    P*TH.     MATRIX 
TERMAT    BECOMES    THE    OUTPUT     MUL T I TERM IN AL    CONNECTIVITY    MAT- 
RIX   A^TER    N-l    COMPUTATIONS. 


MAT{ 35,35) 


C 
C 

C 


c 
c 
c 


c 
c 
c 


SUBcOUTP,c     MATRXT 

IMPLICIT    INTFGFR(A-Z) 

DIMENSION    KMAT(35,35),     FLOMA T ( 3 5 , 35) , 

DIMcNSION    TERMAT  (35  ,3C>) 

DIMENSION    VRTX(35,35)t     NOTVTX ( 3 5 , 35) ,     N0DE(35) 

DIMENSION    NX (35),     NQTNX(35) 

PIMFNSI ON    NLI ST( 35) 

COMMON     /ARR!/     KMAT,FLOMATtCMAT, VPTX, NOTVTX 

COMMnN    /ARR?/     T"MAT 

COMMON    /CHI/    NODE, NLI ST 

COMMON     /TH2/    NX,NOTNX 

COMMON    /PT1/    N,  SHURCS  SINK,  ROW,  MARK,  FLAG,  BRCUT 

COMMON     /PT2/    C0UNT1 ,:0UNT2 ,NCHUNT 

ENTER    CUT    VALUE 

DO    17    1=1, COUNT  1 

DO    16    J=1,C0UNT2 
KUT1=NX ( I ) 
KUT2=N0TNX( J) 
TERMAT (KUT1, KUT 2 )= BRCUT 

TERMAT(KUT2,KUT1 )= BRCUT 

COMPARE    CUT    VALUC    WITH    0THPR    ENTRIES 

DO  15    LL=1,N 

IF  (C0UNT1.EQ.1)     GO    TO    7 

IF  (C0UNT2.E0.1 )     GO    to    9 

IF  (VJODE(LL)  .NE.O)     GO    TO    15 

IF  (FLAG.EQ.O)     GO    TO    8 

CONDENScD    VRTX    IS    ON    THE    SOURCE    SIDE     OF    THE    CUT 

DO    I    KK=1 , NCDUNT 

IF     (NLI ST(KK) .EQ.LL)     GO    TO    2 
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1  CONTINUE 
C 

GC    TO    3 

2  X=NLIST<1) 

IF    (NODE(X) .EQ.O )    GO    TO    8 
GO    TO    15 
C 

c  „  T         CONDENSED    NOTVTX    IS    ON    THE     SOURCE    SIDE    OF    THE    CUT 

3  IKEEP=ROW 

4  IF    (1KEEP.EQ.O)     GO    TO    8 
C 

DC    5    KK=1,N 

IF    (NQTVTX(iKEEP,KK).EQ.LL)     GO    TO    6 

5  CONTINUE 
C 

IKEEP=IKEEP-1 
GC  TO  4 

6  XBAR=NOTVTX( IKEEP*1) 

IF  (NOOE(XBAR).EQ.O)  GO  TO  8 
GC  TO  15 

7  IF  (NODE(LL) .NE.O)  GO  TO  14 

8  IF     (TERNAT(KUT1,LU.LE.TERMAT(KUT1,KUT2)  )  GO  TO  15 
TERMATl KUT1 ,LL) =BRCUT 
TERMAT(LL,KUT1)=BRCUT 

GC  TO  15 

9  IF  (NODE(LL) .NE.O)  GO  TQ  13 
C 

C  DISCARD  THIS  VERTEX  IF  COMPONENT  OF  SCME  CONDENSED 

C  VERTICES 

C 


c 

c 
c 


DC  10  JJ=1,NC0UNT 

IF    (NLI STC J J  J .EQ.LL)     GO    TO    13 

10  CONTINUE 

1KEEP=R0W 

11  IF     (  IKE  EP.  EQ.O)     GO    TO     15 

DC    12    JJ=ltN 

IF     (NOTVTXUKEEP.JJ). EQ.LL)     GO    TO    13 

12  CONTINUE 

IKEEP=IKEEP-1 

IF    (  IKEEP.EQ.O)    GO    TO     15 

GO    TO    11 

13  IF     (TERMAT(KUT2t  LU.LE.TERMAT(KUT1,KUT2)  )     GC    TO    15 
TERMATl KUT2tLL)=BRCUT 

TERMAT(  LUKUT2)  =BRCUT 
GC    TO    15 

14  IF     (C0UNT2.EQ.1)     GO    TO    13 

15  CONTINUE 

16  CONTINUE 

17  CONTINUE 


RETURN 
END 
C 

c 

C     .     . SUBROUTINE    MATRXC 

C 

C     THIS  SUBROUTINE  CONDENSES  A  GRGUP  OF  VERTICES  CN  ONE 
C  SICE  OF  A  CUT  INTO  A  GENERALIZED  SINGLE  VERTEX  AND  OB- 
C  TAINS  A  CAPACITY  MATRIX  CMAT  REPRESENTING  THE  NEW  GRAPH 
C  FORMED. 


C 


SUBROUTINE  MATRXC 

IMPLICIT    INTEGER'(A-Z) 

CIMENSION    KMAT(35f35),    FLOMATt 35 , 35) ,    CM4T(35,35) 

DIMENSION    VRTX(35t35)»     NOT VTX (35 t 35) t     N0DE(35) 

DIMENSION  NLIST(35) 

CCMMON  /ARR1/  KKAT, FLOMAT t CMAT, VRTXt NOTVTX 
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COMMON    /CHI/    NODE,NLIST 

fRIK    4EJK    N, SOURCE, SINK, ROW, MARK, FLAG, BRCUT 
CGMMON    /PT2/    COUNT  1 ,C0UNT2 , NCOU NT 

C  INITIALIZE    CAPACITY    MATRIX    TO    ZERO 


C 

DO    2    1=1, N 
C 

DC    1    J=1,N 
CMAT( I , J)=0 

1  CONTINUE 
C 

2  CONTINUE 
C 

C 

C  FILL    UP    CAPACITY    MATRIX 

C 

DO    13    1=1, N 

IF    (NODE( I ) .EQ.O)     GO    TO    13 

IF     (NCGUNT.EQ.O)     GO    TO    4 
C 

CO    3    LL=l,NCOUNT 

IF     (I .EQ.NLIST(LL) )    GO    TO    13 

3  CONTINUE 

4  NFLAG=0 
JFLAG=0 

DC  12  J=1,N 

IF  (KMAT( I  ,J) .EQ.O)  GO  TO  12 
IF  (NODE(J) .EQ.O)  GO  TO  5 
CMATU  ,  J)~KMAT(  I  ,  J) 
GO  TO  12 
C 

5  INGEX=RCW 

IF    (MARK.EQ.l)    GO    TO    6 
IF     (FLAG.NE.O)    GO    TO    6 
GO    TO    11 
C 

6  CO    7    MM=l,NCOUNT 

IF    (J.EQ.NLIST(MM) i     GO    TO    8 

7  CONTINUE 
C 

GC    TO    9 
C 
C  PROCESSED    VERTICES    (NLIST)     TO    BE    CONDENSED 

8  MFLAG=MFLAG+i 
JTEMP=NLISTC1) 
CKAT( I , JTENP)=MFLAG 
CNAT( JTEMP,I )=MFLAG 
GO    TO    12 

C 

C  UNLABELED  VERTICES  TO  BE  CONCENSED 

C 

9  DO  10  KK=1,N 

IF  (J.EQ.NOTVTX( INDEX, KK) )  GC  to  n 

10  CONTINUE 
C 

INDEX=INDEX-1 

IF  ( INDEX. EQ.O)  GO  TO  12 

GC  TO  9 

11  JFLAG=JFLAG+1 
JTEMP=NOTVTX( INDEX, 1) 
CPATU  ,  JTE^P)=JFLAG 
CNAT( JTEMP,I )=JFLAG 

12  CONTINUE 


C 

c 


CONTINUE 

RETURN 
END 
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C    SUBROUTINE    PKTMAT       ....     

C 

C  THIS    SUBROUTINE    PRINTS    TERMAT    IN    PROPER     FORf    AND     IN 

C    FORMAT    3513. 

C 

C    IF  ARRAY  LENGTH  IS  .IE.  35  A  SINGLE  LOOP  SECTION  PRINTS 

C  THE  MATRIX,  OTHERWISE  IT  PRINTS  THE  FIRST  35  COLUMNS  FOL- 

C  LOkED  BY  THE  NEXT  35  COLUMNS  AND  SO  ON  IN  MULTIPLES  OF  35 

C  USING  MULTIPLE  LOOPS.  THE  REMAINDER  COLUMNS  ARE  PRINTED 

C  IMMEDIATELY  BELOW  THE  LAST  35-CGLUMN  GROUP. 

C 

SUBROUTINE    PRTMAT 

IMPLICIT    INTEGER(A-Z) 

DIMENSION    TERMAT(35,35) 

COMMON    /ARR2/    TERMAT 

COMMON    /PT1/    N, SOURCE, SINK,  ROW,  f'ARK,  FLAG,  BRCUT 
C 
C  HOW    MANY    MULTIPLES    OF    35     IS    N? 

MULPLE=N/35 

IF    (MULPLE.GE  .1)     GO    TO    2 
C 
C  THIS    SECTION    PRINTS    ARRAY    WITH    LENGTH    . LT .    35 


C 


c 
c 


DC    1    I=1,N 

WRITE     (6,8)     I  ,(TERMAT(]  ,J)  , J  =  l ,N) 

CONTINUE 


GO    TO    7 
C 

C  THIS    SECTION    PRINTS    ARRAY    WITH    LENGTH    .GE.    35 

C  IN    GROUPS    OF    35 

2  M=l 
K  =  l 

3  MPLS35=M+35 
C 

DC    4     I=1,N 

WRITE     (6,8)     I  ,(TERMAT( I  ,J)  ,  J=M,MPLS35) 

4  CONTINUE 
C 

M=M+35 

IF    (K.ECMULPLEJ     GO    TO    5 

WRITE     (6,9)     K 

WRITE     (6,10) 

K=K+1 

GC    TO    3 

C  THIS    SECTION    PRINTS    REMAINDER    COLUMNS 

5  REMDR  =  N-MULPLEs-35 

IF     (REMDR.LE.O)     GO    TO    7 
WRITE     (6,11) 
WRITE     (6,12) 
C 

DC    6    L=1,N 

WRITE     (6,8)     L,(TERMAT(L,J) , J=M,N) 

6  CONTINUE 


7  RETURN 

8  FORMAT  CO',     T8 ,     12,     '/',     IX,    3  5  13) 

9  FORMAT  (//TljJ'doNTlNUATloN",    13,    2X:__-0F    TERMAT') 

11  FORMAT  (/./TllT^LAST    SECTION~OF    TERKAT'T 

12  FORMAT  (T14,« '  ' 

END 
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B-4.   PROGRAM  MAXCGN-2A 


C    THIS  PROGRAM  COMPUTES  FOR  THE  MAXIMUM  BRANCH-CONNECT- 

C  IVITY  BETWEEN  ANY  TWO  GIVEN  VERTICES  OF  AN  UNDIRECTED 

C  NETWORK.  THE  NETWORK  GRAPH  IS  REPRESENTED  BY  ITS  CGNNECT- 

C  ION  MATRIX  FOR  COMPUTATION.  THE  NETWORK  CAN  HAVE  PARALLEL 

C  BRANCHES  BETWEEN  VERTICES. 

C 

C    THE  PROGRAM  STARTS  BY  READING- IN  THE  NUMBER  OF  VERTICES, 

C  THE  SOURCE  AND  SINK  VERTICES,  AND  THE  CONNECTION  MATRIX. 

C  THE  CONNECTIVITY  OF  A  CUT(X,X3AR)  IS  DETERMINED  USING  THE 

C  MCCIFIEQ  FORD  AND  FULKERSON  LABELING  ALGORITHM. 

C 

C    THE  OUTPUT  CONSIST  OF  THE  VALUE  OF  THE  CONNECTIVITY  AND 

C  THE  VERTEX  PAIRS  JOINING  EACH  BRANCH  CUT. 

C 

C  THIS    PROGRAM    MAYBE    PREFERRED    OVER    NAXCON-2     WHEN    ONLY 

C    SEVERAL    SOURCE-SINK    PAIRS    ARE    TO    BE     INVESTIGATED. 

C 

C    THE  PROGRAM  DIMENSION  AND  READ  STATEMENTS  MLST  BE 

C  ADJUSTED  FOR  LARGER  NETWORKS. 

C 

C 

IMPLICIT     INTEGER(A-Z) 
C  DIMENSION    KMAT(60,60),     FLOMAT (60 , 60 ) ,     NGCE(60)»     NX(20) 

DIMENSION    N0TNX(2Q),     JSTGRE(20),     ISTCRE(2C) 
C 

READ    (5,17)     N 

READ    (5,18)     (  (KMATU,  J)  ,J  =  1  ,60)  ,I=1,N) 

WRITE  (6,20) 
1  READ  (5,19)  SOURCE, SINK 

IF  (SOURCE.EQ.O)  GO  TO  16 
C 
C  INITIALIZE  FLOW  MATRIX  TO  ZERO 

C 

DC  3  I=1,N 


DO    2    J=1,N 

FLOMAT (  I,  J)=0 

2 

CONTINUE 

C 

3 

CONTINUE 

c 
c 

c 

ERASE     VERTEX    LABELS 

C 

4 

DC    5    J=1,N 
NODE( J) =0 

5 

CONTINUE 

c 

c 

I=SOURCE 

L=0 

M  =  0 

Z  =  C 

BRCUT=0 

NCDE(I) =SOURCE 

C  SEARCH  FOR  AUGMENTING  PATH 
C 

6  DO  7  J=1,N 

C  TEST  FOR  CONNECTING  BRANCH  BETWEEN  VERTEX-PAIR 

IF  (KMAT(I,J).EQ.O)     GO    TO    7 

IF  (FLOMAT (I , J) .  EQ.KMA1  (  I, J  )  )    GO    TO    7 
C 

C  IS    NODE(J)    LABELcD? 

IF  (NODE(J) .NE.O)     GO    TO    7 

NCDE( J) =1 
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C  IS    THIS    THE    SINK    VERTEX? 

IF    (J.EQ.SINK)    GO    TO    10 

JSTORE(M}=J 
7    CONTINUE 
C 

IF    (Z.NE..O)    GO    TO    9 
IF     (M.EQ.0)     GO    TO    11 
C 

C  TRANSFER  CONTENTS  OF  JSTORE  TG  ISTORE 


C 


C 


DC  8  S  =  1,M 

ISTOREi  S)=JSTORE( S) 
8  CONTINUE 


Z  =  M 

M  =  0 
9  I=ISTORE(Z) 

Z=Z-1 

GO  TO  6 
C 

C  WE  REACH  HERE  IF  AN  AUGMENTING  PATH  IS  FOUND 

C  ADJUS7  FLOW  ON  ALL  BRANCHES  ALONG  THIS  PATH. 

10  I=NODE(J) 

FLCMAT(  I,J)=FLCMAT( I , J)+l 

J=I 

IF  (J. EQ. SOURCE)  GO  TO  4 

GO  TO  10 
C 

C  THIS  STEP  IS  REACHED  IF  LABELING  AND  FLOW  AUGMENT- 

C  ATION  IS  COMPLETED.  SEARCH  FCR,  BRANCH  CUTS. 


11  DC  14  I=1,N 

I F  (NODE( IJ.EQ.O)  GO  TO  14 


DC  13  J=1,N 

ADD=1 

IF  (KMAT< I, J ) .EQ.C)  GO  TO  13 

IF  (NGDE( J) .NE.O)  GO  TO  13 
C 
C  IF  FLOW  OCCURS  IN  BOTH  DIRECTIONS,  NET  FLOW  IS  ZERO 

IF  (FLGMATU  ,  J)  .EQ.FLCMAT  (J,  I)  )  GO  TO  13 

IF  (FLOMATU.  ,  J)  .EQ.l)  GO  TO  12 

ACD=FLOMAT(I, J) 

12  8RCUT=BRCUT+ADD 
L  =  L  +  1 

NX(L)=I 
NCTNX(L)=J      • 

13  CCNTINUE 
C 

14  CONTINUE 
C 

c 

C  COMPUTATION    TERMINATES    SINCE    NO    MORE    AUGMENTING 

C  PATH    IS    FOUND. 

IF    (BRCUT.NE.O)     GO    TO    15 

WRITE  (6,21)  SOURCE, SINK, BRCUT 

GC  TO  1 

15  WRITE  (6,22)  SOURCE , SI NK, BRCUT , < NX < J) ,NOT NX ( J ) , J= 1 , L ) 
GO  TO  1 

16  STOP 
C 

17  FCRMAT     (16) 

18  FORMAT     (35I2,/25I2J 

20  FCRMAT     ('0',     10X,      'SOURCE',     10X,      'SINK',     1CX,     «BR- 
1  CONNECTIVITY',     10X,     ' BR-CUTSET') 

21  FORMAT     CO",     12X,     12,     I3X,     12,     17X,    12) 

22  FCRMAT     CO',     12X,     12,     13X,     12,     17X,     12» 
1  (T67,     13,    IX,     '--' ,     13/)) 

END 


77 


B-5.       PRCGRAM    MAXCGN-3 

C  THIS    PROGRAM.    COMPUTES    FOR    THE    MULT  ITERMINAL       VERTEX- 

C    CONNECTIVITY    OF    A    DIRECTED    OR    MIXED    NETWORK.     THE    NETWORK 

C    IS    REPRESENTED    BY     ITS    CONNECTION    MATRIX    FOR    COMPUTATION. 

C    THE    NETWORK    MUST    HAVE    ONLY    SINGLE     BRANCH    BETWEEN    VERTICES 

C 

C    THE  PROGRAM  PROCEEDS  EY  READING-IN  THE  INPUT  CATA  CGN- 

C  SI  STING  OF  THE  TOTAL  NUMBER  OF  VERTICES  AND  THE  CCNNECT- 

C  ION  MATRIX.  IT  THEN  DETERMINES  THE  INITIAL  S'-T  CUT  WITH 

C  THE  FIRST  AND  NTH  VERTICES  AS  SOURCE  AND  SINK  RESP. 

C  HAVING  FOUND  AN  S'-T  CUT,  IT  THEN  FINDS  OTHER  S-T  CUTS 

C  KEEPING  THE  SINK  VERTEX  CONSTANT  AND  VARYING  ONLY  THE 

C  SOURCE  VERTEX  FROM  2  TO  N.  AFTER  COMPLETING  THIS  LOOP, 

C  TFE  SINK  VERTEX  IS  SET  TO  N- 1  AND  THEN  FINDS  AN  S'-T  CUT 

C  BETWEEN  VERTICES  1  AND  N-l .  THE  SOURCE  IS  AGAIN  VARIED 

C  FROM  2  TO  N-l  FOR  THE  OTHER  S-T  CUTS. 

C 

C    THE  PROCESS  IS  REPEATED  UNTIL  N(N-l)  CUTS  ARE  DETER- 

C  MINED.  EACH  REPETITION  INVOLVES  SUBROUTINE  FLCVAR. 

C 

C    THE  OUTPUT  IS  A  PRINTOUT  OF  THE  NXN  VERTEX-CONNECTIVITY 

C  MATRIX. 

C 

C  .......  .  IDENTIFIER  NOMENCLATURE  ... 

C 

C  ARRAYS 

C  KMAT INPUT    CONNECTION    MATRIX    WITH    ENTRIES    0,      It    -1 

C  DMAT S»-T    FLOW    MATRIX    REPRESENTING    AN    S'-T    FLOW 

C  PATTERN.    THIS    MATRIX    IS    UPDATED    BEFORE    A    NEW    SET 

C  OF     S    AND    S«     VERTICES    ARE    DESIGNATED. 

C  FLOMAT-     S-T    FLOW    MATRIX    REPRESENTING    S-T     FLCW    PAT- 

C  TERN    RESULTING    FROM    AUGMENTATION    AND    HOMING    STEPS. 

C  NONZERO    ENTRIES    OF    THIS    MATRIX    UPDATE     CMAT  . 

C  TERMAT-    MATRIX    DENOTING    THE    MAXIMUM    BRANCH-CCNNECT- 

C  IVITY    BETWEEN    ANY    VERTEX-PAIR. 

C 

C  VECTOR 

C  NCDE LIST     OF    VERTICES    OF    THE    NETWORK    GRAPH. 

r- 

C  ADDITIONAL  INFORMATION  ......... 

C 

C  THIS  PROGRAM  IS  CODED  IN  FORTRAN  LANGUAGE.  NETWORK  OF 

C  ANY  SI7E  CAN  BE  HANDLED  BY  THIS  PROGRAM  WITHIN  THE  LIMIT- 

C  ATICNS'OF  THE  COMPUTER  STORAGE.  THE  DIMENSION,  READ,  AND 

C  WRITE  STATEMENTS  MUST  HOWEVER  BE  ADJUSTED  ACCORDINGLY. 

C 

c 

C  MAIN  PROGRAM  .....  

C 

DIMENSION  KMAT(25,25).  FLCMAT( 25  ,25)  ,  DMAT (25, 25) 

DIMENSION  TERMAT ( 25,25) 

DIMENSION  NODE (25) 

CCMMON  /ARR1/  KMA T  ,  FLOMAT , DMAT 

COMMON  /ARR2/  TERMAT 

C CMMUN  /P  T /  SOUR  C  E , SP R I  ME , S I NK , A  MARK , V  T XC  LT  , N , G , F L AG 
C 

READ  (5,11)  N      ,    it    f  „cl  T  ,  fc-» 

READ  (5,12)  UKMATU,  J)  ,J  =  1,25)  ,1=1, N) 

C  INITIALIZE  TERMINAL  MATRIX 


C 


DC  2  X=1,N 

DO  1  Y=l,25 
TERMAT (X,Y)=0 
CONTINUE 
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c 
c 
c 

c 


c 
(. 
c 


c 
c 
c 


c 

c 


c 
c 
c 


c 
c 


c 
c 
c 


2    CONTINUE 

COMPUTE    FOR    THE    N(N-l)     VERTEX    CUTS 

DC    10    K=lfN 
J=N+1-K 

SET    TO    ZERO    THE    FLOW    MATRICES 
DC    4    L=1,N 


DO  3  M=1,N 

CMAT(L,M)=0 

FLCMAT(L,M)=0 

3  CONTINUE 

c 

4  CCNTINUE 

c 

SPRIME=0 

FLAG=0 

DO    9    I=ltN 

IF    (I -EQ.J)    GO    TO    9 

IS    THERE    A    DIRECT    CONNECTION    BETWEEN    SCURCE    AND 
S  INK? 
IF    (KMATU  ,J)  .GT-0  )    GO    TO    5 
SOURCE=I 
SINK=J 
V1XCUT=0 
CALL    FLOVAR 

ENTER    CUT    VALUE    INTO    THE    TERMINAL    MATRIX 
TERMA'U  I  ,J)=VTXCUT 
SPKIME=  I 
GC    TO    6 

5  TERKATl  I  ,JJ=0 
GC    TO    9 

TRANSFER  FLOMAT  ENTRIES  INTC  DMAT 

6  DO  8  L=1,N 

DO  7  M=1,N 

IF  (FLOMAT(LtM).EQ.O)  GO  TO  7 

DNAT(L,M)=FLCMAT(L,M) 

RESET  FLOMAT  TO  ZERO 
FLOMAT (L,M)=0 

7  CCNTINUE 

8  CONTINUE 

FLAG=1 

9  CCNTINUE 

10    CONTINUE 

COMPUTATION  COMPLETED 
WRITE  (6,13) 
WRITE  (6,14) 
CALL  PRTMAT 
STCP 


11  FORMAT 

12  FCRMAT 

13  FCRMAT 

14  FORMAT 
END 


(16) 

M0'^T14,«  VERTEX-CONNECTIVITY    MATRIX*  ) 


(T14, ' 


I 
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C SUBROUTINE  FLOVAR •  •  , 

C 

C    THIS  SUBROUTINE  DETERMINES  AN  S-T  CUT  ONCE  AN  INITIAL 

C  S'-T  CUT  IS  FOUND.  IT  MAXIMIZES  AN  S-T  FLOW  USING  THE 

C  EXISTING  S'-T  FLCW  PATTERN. 

C 

SUBROUT INE    FLOVAR 

IMPLICIT     INTEGER {A-Zi 

DIMENSION    KMAT( 25,25),     FL OM AT ( 25 , 25) t    DMAT(25,25) 

DIMENSICN    NODE(25) 

CCMMON  /ARR1/  KMAT , FLOMAT T DMAT 

CCMMON  /CH/  NODE 

CCMMON  /PT/  SOURCE, SPR  IME ,S  INK, AMARK, VTXCUT ,N ,G,  FLAG 
C 

C  ERASE  ALL  VERTEX  LAbELS 

C 

1  DO  2  J=1.N 
NCDE( J)=0 

2  CONTINUE 
C 

c 

C  SEARCH  FOR  AUGMENTATION  PATH 

CALL  VTXCON 
IF  (AMARK)  3,24,15 

C  AN  S-G  AUGMENTATION  PATH  IS  FOUND 

C  SEARCH  FOR  HOMING  PATH 

3  Q=G 
C 

4  DC  5  J=1TN 

IF  (DMAT(Q,J) .EQ.l)  GO  TO  6 

5  CONTINUE 
C 

GC  TO  1 

6  NODEU)=Q 
G  =  J 

IF  (Q.EQ.GJ  GO  TG  19 

IF  (J.EQ.SINK)  GO  TO  7 

GC  TO  4 
C 

C  A  HOMING  PATH  IS  FOUND 

C  SEARCH  FOR  TRUNCATION  PATH 

7  Z=G 

IF  (G.EQ.SPP IME)  GO  TO  12 

C 

8  DC  9  I=lrN 

IF  (DMAT(  I,Z1  .EQ.l)  GO  TO  10 

9  CONTINUE 
C 

IF  (Z.EO.G)  GO  TO  12 
GC  TO  11 
1C  NODE! I) -I 

IF  (Z-EQ.G)  GC  TO  21 

IF  (Z.EQ.SPRIME)  GO  TO  11 
GO  TO  8 

C  A  TRUNCATION  PATH  IS  FOUND 

C  TRUNCATE  SPRIME-G  PATH 

11  I=NODE(Zi 
DMAKZt  H=0 

IF  (Z.EQ.G)  GO  TO  12 
GO  TO  11 

C  TRACE  AND  ERASE  G-T  PATH 

12  Q=SINK 

13  J=NODE(Q) 
FLCMATt J,QJ=1 
DMAT( J,Q)=0 

Q  =  J 

IF  (J.EQ.G)  GO  TO  14 
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GO  TO  13 
C 
C  AUGMENT  S-G  PATH 

14  IF  (G.EQ. SOURCE)  GO  TO  23 
J=G 

GC  TO  16 
C 
C  AUGMENT  S-T  PATH 

15  J=SINK 

16  I=NODE(J) 

IF  ( I .LT.O)  GO  TO  17 
C 

C      r.«wYI?IEX..IS  WEAKLY  LABELED.  INCREMENT  FWD  FLOW  BY  1 

FLCMATl  I  t  J) =1 

GO    TO    18 
C 
C  VERTEX    IS    STRONGLY    LABELEC    DECREMENT    EACKFLOH    BY    1 

FLOMAKJf  I  )  =  0 

18  J  =  I 
C 

C  IS    SOURCE    VERTEX    REACHED? 

IF    ( J. EQ. SOURCE)     GO    TO    23 

GO    TO     16 
C 
C  A    LOOP    OF    FLOW    IS    FOUND    BY    HOMING    ROUTINE 

19  0  =  G 

20  J=NODE<Q) 
DFAT( J,Q)=0 
Q  =  J 

IF  (Q.EO.G)  GO  TO  1 

GC  TO  2  0 
C 
C  A  LOOP  OF  FLOW  IS  FOUND  BY  TRUNCATION  ROUTINE 

21  Z=G 

22  I=NODE(Z) 
DMATCZy I)=0 
Z=I 

IF  (Z.Et).S)  CO  TO  1 
GC  TO  22 

23  VTXCUT  =  V7XCUTH 
GC  TO  1 

C 

C  NO  MORE  AUGMENTATION  PATH  IS  FOUND.  FLCW  PATTERN 

C  IS  MAX 

24  RETURN 
END 

C 
C 

C  SUBROUTINE  VTXCON   ......... 

C 

C  THIS    SUBROUTINE    LOCATES    AN    S-G    OR    S-T    AUGMENTATION    PATH 

C    USING    FRISCH'S    VERTEX-PAIR    CONNECTIVITY    ALGGRITFM. 

C 

C  IT     IS    INVOKED    BY    SUBROUTINE    FLOVAR    WHICH    SPECIFIES    THE 

C    SOURCE    AND     SINK    VERTICES.     THIS    ROUTINE    LABELS    THE    VERTI- 

C    CES    WEILL    LOCATING    A    PATH.     INFORMATION    WHETHER    THERE     IS 

C    AN    AUGMENTING    PATH    OR    NOT     IS    RETURNED    TO    ROUTINE    FLOVAR. 

C 

SU6R0UT  INE    VTXCON 

IMPLICIT    INTEGER(A-Zi 

DIMENSION    KMAT( 25,25),    FLOMAT( 25  ,  25)  ,    DMAT(25,25) 

DIMENSICN    N0DE(25),    JSTCRE(20),     1ST0RE(20) 

CCMMON     /ARR1/    KMAT  ,  FLCMAT , DMAT 

CCMMON    /CH/    NODE 

CCMMON    /PT/    SOURCE, SPRIME, SINK,  AMARK,VTXCUT,N,G, FLAG 
C 
C  INITIALIZE    VARIABLES 

I = SOURCE 

L  =  0 

M=0 

Z=0 
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c 
c 
c 

c 

c 


c 
c 


c 
c 
c 


c 

c 


c 

c 


c 
c 
c 


c 
c 


c 

c 


c 
c 


c 

c 
c 
c 


AMARK=0 
NCDE(I)=SOURCE 

SEARCH  FOR  AUGMENTING  PATH 

1  DG  6  J  =  ltN 

TEST  FOR  CONNECTING  PATH  BETWEEN  VERTEX-PAIR 
IF  (FLCMATC J,I).EQ.G)  GO  TO  3 
I F  (NQD£( I) .LT.O)  GO  TO  2 

A  BACKWARD  BRANCH  IS  ENCOUNTERED. 
IF  (J  .EC}. SOURCE)  GO  TO  7 
IF  (NODE(J) .LT.O)  GO  TO  6 

VERTEX(J)  IS  EITHER  UNLABELED  OR  WEAKLY  LABELED. 
STRONGLY  LABEL  THIS  VERTEX. 
NCDE(J)  — I 
M  =  M  +  1 

JSTORE(M)=J 
GO  TO  7 

2  IF  (J. EU. SOURCE)  GO  TO  6 
NODE( J) =-I 

GO  TO  5 

3  IF  (KMATU  ,J).NE.l)  GO  TO  6 

A  FORWARD  FLOW  IS  ENCOUNTERED.  TEST  FOR  ZERO  FLOW. 
IF  (FLAG.EU.G)  GO  TO  4 
IF  (DMAT(  I  ,J)  .NF..O)  GO  TO  10 

4  IF  (FLOMATU  ,  J)  .NE.O)  GO    TO  6 
IF  (NODb( Ji.NE.O)  GO  TO  6 
NODE( J  1=1 

IS  SINK  VERTEX  REACHED? 
IF  (J.EQ.SINK)  GO  TO  11 

5  M  =  M+1 
JSTORE(M)=J 

6  CONTINUE 

7  IF  (Z.NE.  0)  GO  TO  9 
IF  (M.EQ.O)  GO  TO  12 

TRANSFER  CONTENTS  CF  JSTGRE  INTO  ISTORE 

DG  0  S=ltM 

1ST0RE(  S)=JSTCR.E  (S) 

8  CUNTINUE 

M  =  0 

9  I=ISTORE(Z) 
Z  =  Z-1 

GC  TO  1 

AN  S-G  AUGMENTATION  PATH  IS  FCUNO 

10  G  =  I 

AMARK=-1 
GC  TO  13 

AN  S-T  AUGMENTATION  PATH  IS  FOUND 

11  ANARK=1 
GO  TO  13 

NO  AUGMENTATION  PATH  FOUND 

12  AMARK=0 

13  RETURN 
END 
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£     SUBROUTINE    PRTMAT       

£    CrruiyiL?^ER0UTINE    POINTS    TERMAT    IN    P-ROPER    FCRN    AND     IN 
L     hi. K MAT    35  13  • 

C 

r  tuJFma?2tv     Lri*§IB.,IT§,-'^*     35    A    SINGLE    LOOP     SECTION    PRINTS 

r  ^r1Ern^J8I^cCIKfWiIE     IT     PRINT^    THE     FIRST    35     COLUMNS    FOL- 
r  S,X,  1^1,  £E,XT    ^5    COLUMNS    AND    SO    ON    IN    MULTIPLES    OF    35 

9  t  §ic^T5JykriPti,  k^'°?S?     THE    REMAINDER    COLUMNS    ARE    PRINTED 

C  IMMEDIATELY    BELOW    THE    LAST    35-COLUMN    GROUP. 

SUEROUT INE    PRTMAT 

IMPLICIT    INTEGER(A-Z) 

DIMENSION    TERMAT(25,25) 

CCMMON    /ARR2/    TERMAT 

COMMON    /PI/    SOURCE, SPR IME , SI NK, AMARK , VTXCUT  ,N ,G , FLAG 

C  HOW    MANY    MULTIPLES    CF    35    IS    N? 

MULPLE=N/35 

IF     (MULPLE.GE .1)     GO    TO    2 
C 
C  THIS    SECTION    PRINTS    ARRAY    WITH    LENGTH    . LT.    35 

DC    1    1=1,  N 

WRITE     (6,8)     I, (TERMAT( I ,J) , J=1,N) 

1  CONTINUE 
C 

GO    TO    7 
C 

C  THIS    SECTION    PRINTS    ARRAY    WITH    LENGTH    .GE.    35 

C  IN    GROUPS    OF    35 

2  M  =  l 
K  =  l 

3  MPLS35=M+35 
C 

DC    4    I=1,N 

WRITE     (6,8)     I ,( TERMAT ( I ,J) , J=M,MPLS35) 

4  CONTINUE 
C 

M=M+35 

IF     (K.EQ.MULPLE)     GO    TO    5 

WRITE     (6,9)     K 

WRITE     (6,10) 

K  =  K+1 

GO    TO    3 
C 
C  THIS    SECTION    PRINTS    REMAINDER    COLUMNS 

5  REMDR=N-MULPLE-'-3  5 

IF  (REMCR.LE.O)  GO  TO  7 
WRITE  (6,11) 
WRITE  (6,12) 
C 

DC  6  L=1,N 

WRITE  (6,8)  L,(TERMAT(L,J) , J=M,N) 

6  CONTINUE 

7  RETURN 


C 

c 


8  FORMAT  ('0',  T8,  12,  «/»,  IX,  3513) 

9  FORMAT  (//T14, "CONTINUATION « ,  12,  2X,  'OF  TERMAT') 
1C  FORMAT  (T14,' — ^ — - -•) 

11  FORMAT  (//T14,'LAST  SECTION  OF  TERMAT') 

12  FORMAT  (T14,« •') 

END 
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B-6.        PRCGRAM    MAXC0N-3A 


C  THIS    PROGRAM    COMPUTES    FOR    THE    VERTEX-CONNECTIVITY    OF    A 

C    DIRECTED,    AN    UNDIRECTED,    OR    MIXED    NETWORK    USING    FRISCH'S 

C    VERTEX-PAIR    CONNECTIVITY    ALGORITHM. 

C 

C    THE  PROGRAM  STARTS  BY  READING-IN  THE  INPUT  CATA  CONSIS- 

C  TING  OF  THE  NUMBER  OF  VERTICES,  THE  CONNECTION  MATRIX, 

C  AND  THE  SOURCE  AND  SINK  VERTICES  AS  CHOSEN.  TFENCE  IT 

C  FINDS  AN  AUGMENTING  PATH  USING  THE  MODIFIED  LABELING  AL- 

C  GCRITHM  AND  I NCRE AS ES/ OECR E AS ES  THE  FLOW  ONCE  SUCH  PATH 

C  IS  FUUMD.  THE  PROCESS  IS  REPEATED  UNTIL  NO  MOPE  AUGMENT- 

C  ING  PATH  IS  FOUND.  IT  THEN  SEARCHES  FOR  THE  VERTEX  CUTSET 

C  BY  FINDING  A  LABELED  VERTEX  AND  AN  UNLABELED  CNE  WHICH 

C  ARE  BOTH  CONNECTED  AND  HAVE  UNITY  F LCW  BETWEEN  THEM. 

C 

C    THE  OUTPUT  IS  A  PRINTCUT  OF  THE  SOURCE  VERTEX,  THE  SINK 

C  VERTEX,  AND  A  LIST  OF  VERTICES  COMPRISING  THE  CUTSET. 

C 

C    FCR  NETWORKS  LARGER  THAN  THE  DECLARED  DIMENSION,  THE 

C  DIMENSION  AND  READ  STATEMENTS  MUST  BE  ADJUSTED. 

C 

C 

IMPLICIT  INTEGER(A-Z) 

DIMENSION    KMAT(50,50),    FLOMAT( 50 , 50) 

DIMENSION    NODE150),    JSTORE(20),     ISTORE(20),    VTXCUT(IO) 

READ    (5,24)    N 

READ    (5,25)     ( (KMAT ( I , J ) , J= 1 , 50)  , I  =  1, N) 

WRITE     (6,27) 

1  READ    (5,26)     SOURCE, SINK 
IF     (SOURCE. EQ.O)     GO    TO    23 

r 

C  IS  THERE  A  DIRECT  CONNECTION  BETWEEN  SOURCE 

C  AND  SINK? 

IF  (KMAKSOURCE  ,SINK)  .GT.O)  GO  TO  21 

C 

C  INITIALIZE    FLOW    MATRIX    TU    ZERO 

C 

DO    3    I=lfN 
C 

DO    2    J=1,N 

FLOMATt  I, J  )  =  0 

2  CCN'TINUE 
C 

3  CONTINUE 
C 

c 

C  ERASE    VERTEX    LABELS 

C 

4  DO    5    J=1,N 
NCDE(J)=0 

5  CONTINUE 
C 

C  INITIALIZE    VARIABLES 

I = SOURCE 

L  =  0 

M=0 

Z=0 

NODE( I)=SOURCE 
r 
C  SEARCH    FOR    AUGMENTING    PATH 

C 

6  DO    10    J=ltN 

C  TEST    FOR    CONNECTING    PATH    BETWEEN    VERTEX-PAIR 

IF    (FLOMAKJi  D.tQ.O)     GO    TO    5 
IF    (NODE(I).LT.O)     GO    TO    7 

C  A    BACKWARD    bRANCH    IS    ENCOUNTERED. 
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IF    (J.EQ. SUURCE}     GO    TO    11 
IF     (NODHJ  ).LT.O)     GO    TO    10 

£  VE?I^,(^J     IS    EITHER    UNLABELED    OR    WEAKLY    LABELED. 

C  kinn-(     STRONGLY    LABEL    THIS    VERTEX. 

NODt( J) =-I 

M=M+1 

JSTORE(M)=J 

GO    TO    11 

7  IF    (J.EQ. SOURCE)     GO    TO    10 
NODE( J)=-I 

GO    TO    9 
C 

8  IF     (KMAT(-J,J).NE-1)    GO    TO    10 
C 

C  A  FORWARD  BRANCH  IS  ENCOUNTERED.  TEST  FCR  0  FLCW. 

IF  (FLOMATU,  J)  .NE.O)  GO  TO  10. 

IF  <NQDE( J).NE.O)  GO  TO  10 

NODE(J)=I 
C 
C  IS  SINK  VERTEX  REACHED? 

IF  (J.EQ. SINK)  GO  TO  14 

9  M  =  M+1 

J  STORE ( M) =J 

10  CONTINUE 
C 

11  IF  (Z.NE.O)  GO  TO  13 
IF  (M.EQ.Q)  GO  TO  17 

C 

C  TRANSFER  CONTENTS  OF  JSTORE  INTO  1ST0RE 

C 

DC  12.  S  =  l,M 

ISTORE( S)=JSTORE(S) 

12  CONTINUE 
C 

Z  =  M 
M  =  0 

13  I=ISTORE(Z) 
Z  =  Z-1 

GO  TO  b 
C 

C  WE  REACH  HERE  IF  AN  AUGMENTING  PATH  IS  FOUND. 

C  ADJUST  ALL  FLOWS  mLONG  THIS  PATH. 

C  TEST  LABEL  ON  VERTEX(J) 

14  I=NGDL-(J) 

IF    (I.LT.O)    GO    TO     15 
C 
C  VERTEX    IS    WEAKLY    LABELED.     INCREMENT     FWC    FLOW    BY    1 

FLCMATC I, J)=l 

GO    TO    16 
C 
C  VERTEX     IS    STRONGLY    LABELED.     DECREMENT    BACKFLOW 

15  I=-I 
FLOMAT( J, I )=0 

16  J=I 
C 

C  IS  SOURCE  VERTEX  REACHED? 

IF  (J.EQ. SOURCE)  GO  TO  4 

GC  TQ  14 
C 
C 

C  VERTEX  LABELING  AND  FLOW  AUGMENTATION  IS  COMPLTD. 

C  THIS  SECTION  SEARCHES  FCR  THE  VERTEX  CUTS 

C 

17  DC  20  1=1 tN 

IF  (NODE( I) .EQ.O)  GO  TO  20 
C 

DC  19  J=1,N 

IF  (NODt( J). NE.O)  GO  TO  19 

IF  (FLOMATU  ,  J)  .EQ.O)  GO  TO  19 

IF  (I .EQ. SOURCE)  GO  TO  18 

L  =  L  +  1 
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VTXCUT(  Li  =1 
GO    TO    20 

18  L=L+1 
VTXCUK  L)=J 

19  CONTINUE 
C 

20  CONTINUE 
C 

c 

c 

C  COMPUTATION    TERMINATES    SINCE    NO    MORE    AUGMENTING 

C  PATH    IS    FOUND.    PRINT    CUTSET. 

IF    (L.EQ.O)    GO    TO    21 

WRITE  (6,28)  SOURCE, SINK, (VTXCUT(J)  ,J=1  ,L) 

GC  TO  1 

21  WRITE  (6,29)  SOURCE, SINK 
GC  TO  1 

22  WRITE  (6,30)  SOURCE, SINK 
GO  TO  1 

23  STOP 
C 

24  FORMAT  (16) 

25  FCRMAT  (2513/2513) 

26  FCRMAT  (216) 

27  FORMAT     COS     10X,     'SOURCE',     10X,     'SINK1,     10X, 
1  "VERTEX-CUTSET1) 

28  FCRMAT  CO1,  12X,  12,  13X,  12,  (T42,  413)/) 

29  FORMAT  C 0 5 , 12X, I  2 ,  13X, I  2 , T4 1,  8 NO  FEASIBLE  FLCW1) 

30  FCRMAT  CO',  12X,  12,  13X,  12,  T44,  «08) 
END 
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B-7.   PRCGRAM  MAXC0NK4 


C    THIS  PROGRAM  VERIFIES  THE  R-CONNEC TI VI TY  OF  AN  UNDIR- 

C  ECTED  NETWORK  FOR  A  GIVEN  UNIFORM  REQUIREMENT  R. 

C 

C    INPUT  DATA  ARE  THE  NUM&ER  OF  VERTICES,  THE  REQUIREMENT 

C  R,  AND  THE  NETWORK  CONNECTION  MATRIX. 

C 

C  THE    PROGRAM    READS- IN    THE     INPUT    DATA    THEN    DETERMINES    THE 

C     NUMBER    OF    VERTEX    DISJOINT    PATHS,     STARTING    WITH    VERTEX    1 

C    AS    SOURCE,     TC     EACH    OF     THE    REMAINING    VERTICES.     IT    EMPLOYS 

C    FRISCH'S    ALGORITHM    FOR    VERTEX-PAIR    CONNECTIVITY    AND 

C    KLEITMAN'S    THEOREMS    FOR    OBTAINING    THE    OVERALL    CONNECT  I- 

C    V  ITY. 

C 

C  IF     A    CONNECTIVITY    LESS    THAN    THE    •REQUIREMENT     IS    ENCOUNT- 

C    ERED,     THE    PROGRAM    CEASES    FROM    FURTHER    VERIFICATION    AND 

C    PRINTS    OUT    THE    CONNECTIVITY    VALUE    SPECIFYING    THE    CORRESP- 

C    GNDIN'G    SOURCE    VERTEX.     OTHERWISE    THE     PROGRAM    WILL     INDICATE 

C    THAT    THE    NETWORK     IS    AT    LEAST    R-CONNECTED. 

C 

C 

IMPLICIT  INTEGER (A-Z) 

DIMENSION  KMAT(50,50),  ILIST(IO) 

COMMON     /ARR/    KMAT 

COMMON    /CHAIN/     IL 1ST 

COMMON    /PT/    SOURCE , SI NK,LL, C,N, RTEMP 
C 

READ    (5,7  1    N,R 

READ    (5,8)     ( IKMAT( I  ,J)  , J=l  ,50)  *  1^1  ,N) 

L=0 

INDEX=0 

RTEMP=R 
C 

DC    4    1=1, R 

IPLUS1^I+1 
C 

DC    3    J=IPLUS1,N 

C  LIST    DOWN    PREVIOUSLY    DESIGNATED    SINK     VERTICES    THAT 

C  ARE    CONNECTED    TO    VERTEX    J. 

LL  =  0 

JLESS1=J-1 

IF    ( JLESS1.LT. IPLUS1)     GO    TO    2 

C 

DO    1     L=IPLUS1,JLESS1 

IF    (KMAT(L,J) .EQ.O)     GO    TO    1 

LL=LL+1 

ILIST {LL)=L 

1  CONTINUE 
C 

IF     (LL.GE-RTEMP)     GO    TO    3 
r 
C  SEARCH    FOR    DISJOINT    PATHS 

2  SOURCE=I 
SINK=J 

Q  =  0 

CALL    VTXCON 

S=Q+LL+ INDEX 

IF  (S.GE.R)  GO  TO  3 

GO  TO  5 

C  ARE  ALL  (N-I)  VERTICES  PROCESSED? 

3  CONTINUE 
C 

RTEMP=RTEMP-I 
INDEX=INDEX+1 

4  CONTINUE 
C 

C 
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C  VERIFICATION  OF  RTH  VERTEX  AS  SOURCE  COMPLETED. 

WRITE  (0,9)  R 

GO  TO  6 
C 
C  CONNECTIVITY  IS  LESS  THAN  R.  TERMINATE. 

5  WRITE  ( 6, 10)  S 
WRITE  (6,11)  I 

6  STOP 
C 

7  FORMAT  (2  16) 

8  FORMAT  (25  13/2513) 

9  FORMAT  ( '0' ,T10,  'NETWORK  IS  AT  LEAST*  , I  3 , "-CONNECTED1  ) 

10  FORMAT  («0« ,T10, 'NETWORK  IS  ONL Y « , 13, '-CGNNEC TED  WITH' 

11  FORMAT  CO',  T13,  «  VERTEX4  ,  I  3  ,2  X  ,  »  AS  SOURCE') 
END 

C 
C 

C «  SUBROUTINE  VTXCON  .  .  .  .  - . 

C 

C  THIS    SUBROUTINE    DETERMINES    THE    NUMBER    OF     VERTEX    Dl  S- 

C    JOINT    PATHS     BETWEEN    TWO    GIVFN    VERTICES    USING    FRISCH'S    AL- 

C    GOPITHM.     IT    IS    HOWEVER    MODIFIED    BY     INCORPORATING    KLEIT- 

C    MAN'S    THEOREMS    TO    REDUCE    COMPUTATIONAL    TIME. 

C 

SUBROUTINE  VTXCON 

IMPLICIT  INTEGER(A-Z) 

DIMENSION  KMAT(50,50i,  FLOM AT ( 50  ,  50) 

DIMENSION  N0DE(50),  JST0RE(20),  ISTOREC20),  ILIST(IC) 

CCMMON  /ARR/  KMAT 

COMMON  /CHAIN/  I  LI  ST 

COMMON  /PT/  SOURCE,SINK,LL,Q,N,RTEMP 
C 

C  INITIALIZE    FLOW    MATRIX    TO    ZERO 

C 

DO    2    I=SOURCE,N 
C 

DO    1    J=SOURCE,N 

FLOMATt I , J)=0 
1    CONTINUE 


C 

c 


2    CONTINUE 
MARK=0 


C 

C  ERASE    VERTEX    LABELS 

C 

3  DC    4    J=SOURCE,N 
NODE(J)=0 

4  CONTINUE 
C 

C 

C  INITIALIZE    VARIABLES 

I=SOURCE 

L  =  0 

M=0 

Z=0 

NODE( I )=SOURCE 

SPLUS1=S0URCE+1 

C  IS    SOURCE    CONNECTED    TO    SINK? 

IF     (KMAT (SOURCE, S INK)  .hQ.O)     GO     TO 
KMAT( SOURCE ,SINK) =0 
KMAT (SINK, SOURCE )=C 
Q  =  Q+1 
MARK=1 

C  SEARCH    FOR    AUGMENTING    PATH 

C 

5  DO  11  J=SPLUS1,N 

REMOVE 

C  S  INK 
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C  REMOVE  PROCESSED  VERTICES  THAT  ARE  CONNECTED  TO 


IF    (LL.EQ.O)    GO    TO    7 
C 

DC    6    K=i,LL 

IF    (J.EO. ILIST(K) )     GO    TO    11 

6  CONTINUE 
C 

c 

C  -7     rr     ,IF£L£9R.    CONNECTING    PATH    BETWEEN    VERTEX-PAIR 

7  IF     (FLOMAT ( J, T) . EQ.O)     GO    TO    9 
I  F    (NODE(  D.LT.OJ     GO    TO    8 

C  A    BACKWARD    BRANCH    IS    ENCOUNTERED 

IF     (NODE( J) .LT.O)     GO    TO    11 
IF    (J. EQ. SOURCE)     GO    TO     12 

C  VERTEX(J)     IS    EITHER    UNLA8ELEC    OR    WEAKLY    LABELED. 

C  STRONGLY    LABEL    THIS    VERTEX 

NODE(J)=-I 

M  =  M+1 

JSTORE(M)=J 

GC  TO  12 

8  IF  (J. EQ. SOURCE)  GO  TO  11 
NQDE(J)=-I 

GC  TO  10 

9  IF  (KMAT( I ,J) .NE.1J  GO  TO  11 
C 

C  A  FORWARD  BRANCH  IS  ENCOUNTERED.  TEST  FOR  ZERO  FLOW 

IF  (FLOMATU  ,J).NE.O)  GO  TO  11 

IF  (NODEU ) .NE.O)  GO  TO  11 

NODE* J) =1 
C 
C  IS  SINK  VERTEX  REACHED? 

IF  (J. EU. SINK)  GO  TO  15 

10  M=M  +  1 
JSTORE(M)=J 

11  CONTINUE 
C 

12  IF  (Z.NE.O)  GO  TO  14 
IF  (M.EQ.O)  GO  TO  19 

C 

C  TRANSFER  CONTENTS  OF  JSTORE  TO  I  STORE 

C 

DO  13  S-1,M 

ISTORE(S)=JSTCRE(S) 

13  CONTINUE 
C 

Z=M 
M  =  0 

14  I=--ISTORE(Z) 
Z=Z-1 

GO  TO  5 
C 

C  WE  REACH  HERE  IF  AN  AUGMENTING  PATH  IS  FOUND 

C  ADJUST  ALL  FLOWS  ALONG  THIS  PATH. 

15  Q=Q+1 
U=G+LL 

IF  (U.GE.RTEMP)  GO  TO  19 
C  TEST  LABEL  ON  VERTEX  J 

16  I ---NODE  (J) 

IF  (I.LT.O)  GO  TO  17 
C 
C      VERTEX  IS  WEAKLY  LABELED.  INCREMENT  FORWARD  FLOW  BY  1 

FLOMAT  (  I,  JJ=1 

GO  TO  18 
C 
C  VERTEX  IS  STRONGLY  LABELED.  DECREMENT  BACKFLOW  BY  1 

17  I=-I 
FLCMAT-(  J,  I  )=0 

18  J  =  I 
C 

C  IS  SOURCE  VERTEX  REACHED? 

IF  (J. EQ. SOURCE)  GO  TO  3 
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c 

c 


19 


20 


GC  TO  16 

VERTEX  LABELING  AND  FLOW  AUGMENTATION 
IF  (MARK.EQ.O)  GO  TO  20 
KMAT( SOURCE iS INK) =1 
KMAT( SINK, SOURCE )  =  1 
RETURN 
END 


IS  COMPLETED 
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